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Abstract 

In many-body fermionic systems at finite density correlation functions of the density operator 
exhibit Friedel oscillations at a wavevector that is twice the Fermi momentum. We demonstrate 
the existence of such Friedel oscillations in a 3d gravity dual to a compressible finite-density state 
in a (1+1) dimensional field theory The bulk dynamics is provided by a Maxwell U{1) gauge 
theory and all the charge is behind a bulk horizon. The bulk gauge theory is compact and so 
there exist magnetic monopole tunneling events. We compute the effect of these monopoles on 
holographic density-density correlation functions and demonstrate that they cause Friedel oscilla- 
tions at a wavevector that directly counts the charge behind the bulk horizon. If the magnetic 
monopoles are taken to saturate the bulk Dirac quantization condition then the observed Fermi 
momentum exactly agrees with that predicted by Luttinger's theorem, suggesting some Fermi sur- 
face structure associated with the charged horizon. The mechanism is generic and will apply to 
any charged horizon in three dimensions. Along the way we clarify some aspects of the holographic 
interpretation of Maxwell electromagnetism in three bulk dimensions and show that perturbations 
about the charged BTZ black hole exhibit a hydrodynamic sound mode at low temperature. 
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I. INTRODUCTION 

This paper will be concerned with a holographic description of a particular compressible 
phase of matter. Recall that a compressible state is one in which the density of a U{1) 
charge is a continuous function of various parameters such as the chemical potential /x or the 
temperature T. In traditional perturbative quantum field theory, it is generally understood 
that for such a state to persist to arbitrarily low temperatures without breaking the U{1) 
symmetry the charge cannot be carried by bosons, which have a tendency to condense and 
form a superfiuid phase. The charge then must then be carried by fermionic excitations, 
which will fill out one or more Fermi surfaces. One would then naively expect that in any 
compressible state at sufficiently low temperature various correlation functions - such as 
those of the charge density itself - should exhibit sharp singularities in momentum space 
associated with these Fermi surfaces. 

In fact the locations of these singularities are tightly constrained by Luttinger's theorem 
[1], which relates the volume enclosed by the Fermi surface to the charge density p measured 
in units of the fundamental charge quantum q^.. For example, in a 1 + 1 dimensional field 
theory, this relation reads 

P=|^(2M, (1-1) 

with k = ±kF the locations of the Fermi points. In the standard field theoretical examples 
the density-density correlator {p{k)p{—k)) exhibits singularities at A; = 2kF- Upon Fourier 
transformation the resulting cos(2/ci?x) structures in position space are called Friedel oscil- 
lations. We see that these oscillations are actually rather important: they provide a direct 
probe of the underlying Fermi surface. We note at this point that in (1+1) dimensions 
the above discussion degenerates slightly in that there is really no fundamental difference 
between a density of bosons and of fermions: there is no true superffuid phase in one spa- 
tial dimension and a finite density of bosons also exhibits Friedel oscillations at the same 
wavevector (1.1). Thus the correlation between Friedel oscillations and charge density is 
even tighter in (1+1) dimensions. 

A. Problem: the nature of horizon charge 

It is thus somewhat perplexing that gauge/gravity duality [2] presents us with many 
examples of compressible states which at first glance do not exhibit any such singularities. 
Here we study a strongly coupled field theory with a large number ('W^") of degrees of 
freedom. The dual gravitational description of any such state involves a charged black hole 
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horizon that sources electric flux in the bulk. This flux then penetrates through the bulk 
to the conformal boundary of the spacetime, where its boundary value can be interpreted 
as the charge density p of the dual fleld theory via the usual AdS/CFT dictionary. In 
some examples this horizon structure persists to zero temperature. Despite extensive study 
of these systems in applications of holography to condensed matter physics (see e.g. [3- 
5] for reviews), it is generally diflicult to identify exactly what fleld-theoretical degrees 
of freedom are carrying the charge. The current understanding of such a system is that 
the charge is carried by gauge-charged excitations in the dual fleld theory [6-10]: in the 
holographic literature such phases are often called "fractionalized" ^ . Note that in such a 
system we cannot easily access any correlators corresponding to fundamental fermions, which 
presumably are not gauge-invariant. However we can compute the density-density correlator, 
and despite extensive study so far no correlators computed in these backgrounds demonstrate 
any special structure in momentum space [11-13]^. Clearly there is some tension with the 
fleld-theoretical intuition described above. 

It is the goal of this paper to address this puzzle in the context of a speciflc model, of the 
3d gravitational dual of a 1+1-dimensional quantum liquid. We will demonstrate that if we 
take into account appropriate non-perturbative effects in the bulk, the expected singularities 
are indeed present in the density-density correlation at the expected location in momentum 
space, although their amplitude is strongly suppressed. This suggests - as one would have 
hoped ^ that some dual Fermi surface structure exists in these charged black holes, despite 
the fact that all of our calculations are of gauge fleld dynamics on a 3d curved background, 
with no explicit charged fermions in sight. 

We stress that this problem is different from the inclusion of explicit fermions in the bulk. 
Such bulk fermions are dual to gauge-invariant charged operators in the fleld theory. These 
systems have been extensively studied and often possess explicit Fermi surfaces ([15-19]; see 
[20] for a review) in the correlation functions of the dual fermionic operators. These gauge- 
neutral fermions contribute a parametrically small (0(1) vs. 0{N'^)) charge density that is 
outside the black hole horizon. This is dual to the fact that the charge density is carried by 
gauge-invariant excitations in the fleld theory and so can essentially be understood in the 
framework of traditional fleld theory, as has been emphasized recently by various authors 
[9, 21, 22]. In particular, Friedel oscillations sourced by these explicit bulk fermions have 
been studied in [23]. There are no such bulk fermions in our description. 



^ The name arises from an analogy with related systems in condensed matter theory, where any gauge group 
is emergent and associated with the dissociation of the electron into degrees of freedom with fractional 
quantum numbers. These fractional excitations should be identified with gauge-charged degrees of freedom 
in holography. 

^ Note that singularities in momentum space in a holographic model have been found in the recent work 
[14]; the setup is rather different, involving a density of strings rather than particle-like excitations, and 
the connection to the results of this paper is not clear. 
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B. Monopoles and Berry phases 



We now briefly describe the relevant structures in the bulk, leaving a detailed description 
to the next section. We will study a strongly coupled 1+1 dimensional fleld theory with a 
conserved current j^. This fleld theory has a three-dimensional gravitational dual which we 
take to be weakly coupled. We will assume that the current is dual to a bulk gauge fleld 
Aj^^ and that the leading dynamics of this gauge fleld are given by a bulk Maxwell term 
{dAy. We note that while this is standard in higher dimensional examples of gauge/gravity 
duality generally in AdS3/CFT2 one expects a different structure involving bulk Chern- 
Simons terms. This is not quite the model we will study and we will elaborate on the 
relation between this model and ours in Section II. 

Turning on a chemical potential for the boundary theory current, we flnd that the relevant 
bulk gravity solution is the charged BTZ black hole [24, 25]. The detailed structure of this 
solution is not important. The key fact is that there is a nonzero bulk electric fleld present, 
sourced by the horizon and pointing outwards in the emergent holographic direction: 

F'\r) oc p . (1.2) 

The standard AdS/CFT prescription relates fluctuations of the bulk gauge fleld Aj^ about 
this background to the correlation functions of the boundary theory current j^. As men- 
tioned above, correlation functions previously calculated using this method do not reveal 
any nontrivial structure in momentum space. This is perhaps not surprising. If we take the 
relation (1.1) seriously, we see that the location of any such singularity - i.e. the value of 
2k F - depends not only on p but also on ^g, the charge of a single quantum excitation in 
the fleld theory. The black hole horizon and linearized fluctuations around it know nothing 
of any ^g, and so it seems clear that they cannot reproduce this answer. 

However, does have a bulk interpretation. It is expected that in any theory of quantum 
gravity - such as the one that we are studying in the bulk - all gauge symmetries should be 
compact. In particular Am should be a compact U{1) gauge fleld with a minimum quantum 
of charge, and it is easy to see that this charge quantum should be identifled with the 
quantum of charge characterizing the boundary theory Hilbert space. Thus the deflnition 
of the bulk theory does contain a q^^ even if the black hole solution does not seem to make 
explicit use of it. 

Once we accept that the bulk gauge theory is compact, however, a new ingredient presents 
itself: magnetic monopoles. In three Euclidean dimensions these objects are localized in both 
space and time, and so should be thought of as instantons. If a monopole with magnetic 
charge q^ is placed at a point then we flnd that magnetic flux is created there, i.e. 

dF{X) = qm5^''\X - Xm) (1.3) 

On general grounds we expect these monopoles to represent allowed tunneling events in the 
theory. They can have a profound effect on infrared physics; for example in flat space in 
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three dimensions they drive confinement of compact U{1) gauge theory [26]. In the remain- 
der of this paper we wiU compute the eflFect of bulk monopoles through AdS/CFT on the 
density-density correlation function of the dual field theory. The calculation is conceptually 
quite simple: monopoles in the bulk source the bulk gauge fields, resulting in a nontrivial 
contribution to the boundary theory density-density correlators. As we explicitly demon- 
strate, the problem can be formulated in terms of a Witten diagram as in Figure 1, where 
the monopoles are integrated all throughout the bulk with the appropriate action cost. We 
will see that these monopoles provide a nontrivial probe of the charge density. 



pix) X ► p(0) 




E \ \ \ 

Th 

FIG. 1: Witten diagram used in calculation of contribution of monopole-anti-monopole pair to 
boundary theory density-density correlator. Solid wiggly lines indicate bulk-to-boundary correla- 
tors, and dotted gray lines indicate monopole-antimonopole interaction energy. 

In particular, the boundary theory charge density manifests itself as the bulk background 
electric field (1.2). In the presence of such an electric field, each monopole event is weighted 
by a nontrivial phase^ which we will call a Berry phase in a slight abuse of notation. As we 
will explain in detail later, this phase takes the form 

SsiXm) = iqmpXm (1-4) 

where Xm is the spatial coordinate of the monopole in the direction parallel to the hori- 
zon. This is completely analogous to the Aharonov-Bohm phase accumulated by an electric 
charge in a background magnetic field. We see that the phase associated with the monopole 
oscillates in space with period pg^, indicating a special feature at a particular position in mo- 
mentum space. This feature turns out to be the expected Friedel oscillation. The existence 
of this phase is quite easy to understand: the bulk of the paper is devoted to determining 
the form of the interaction energy between monopoles, which is crucial in determining the 
precise form of the singularity. Similar Berry phases associated with monopole events play 
an important role in traditional condensed matter physics (see e.g. [27, 28]). 
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Comparing to the field-theoretical expression (1.1) we see that there is a nontrivial inter- 
play between electric and magnetic charges. Fascinatingly, we will find precise agreement 
between the field-theoretical and gravity expressions if the spectrum of magnetic and electric 
charges in the bulk saturates the Dirac quantization condition, i.e. qeqm = 27r. We comment 
on this (and associated issues) in the discussion section. 

We now briefiy summarize the paper. In Section II we define the model that we will 
be studying - Maxwell electromagnetism on AdSs and discuss some tree-level aspects: in 
particular we explain how it can be understood as a two-dimensional field theory with a 
conserved current and a logarithmically marginal deformation. In Section III we develop 
a formalism for computing monopole instanton corrections to holographic correlators. In 
Section IV we compute this correction and display the Friedel oscillations in the high- 
temperature limit, in which gravity decouples from the gauge field, greatly simplifying the 
calculation. 

In Section V we turn to the low temperature limit, in which the mixing with gravity is 
crucial. Here we first study the system at tree level and exhibit a sound mode in fiuctuations 
about the charged BTZ black hole. In addition to being interesting in its own right, this 
sound mode determines the long-range correlations between monopoles and so is needed 
for the instanton calculation. We then present two routes to the zero-temperature Friedel 
oscillations: Section VB is very quick, but we feel that it obscures some important physics 
involving the interaction of monopoles with gravity, which we develop in detail in Section 
VC. We summarize our results and discuss some implications in Section VI. The trusting 
and impatient reader can skip to the conclusion section with little loss of continuity. 



II. PRELIMINARIES 

In this section we define our model and establish some facts that will be needed to derive 
the Friedel oscillations. We will study a 1+1 dimensional field theory with a conserved 
current j^. We assume that this theory has a weakly coupled gravity dual governed by the 
classical Euclidean action 

where is the bulk gauge field dual to the boundary theory current and F = dA. /x, u run 
over the three bulk directions, and a, b will run over boundary theory spacetime coordinates. 
In general we will use y to refer to boundary theory coordinates, so = (r, y^) = (r, r, x)^. 
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Our conventions for the Levi-Civita tensor in 2d and 3d are 



Two-index epsilon always lives in flat 2d space and refers to the symbol, not the tensor. 
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We will work with spacetimes which are asymptotically AdS, taking the form 

ds^ = {dr^ + dx^) + ^ (2.3) 

as r ^ oc. The AdS radius has been set to 1 for notational convenience. We will generally 
work in Euclidean signature but will sometimes Wick rotate to real time via 

t —IT uj iuje iS —Se . (2.4) 

In gauge/gravity duality in higher dimension an action such as (2.1) is standard for the 
bulk dual of a field theory with a conserved current. However in two dimensions the story is 
usually different. A conserved current in a two-dimensional conformal theory is associated 
with a higher symmetry: Ward identities guarantee the separate conservation of both the 
holomorphic and antiholomorphic part of the current j^, resulting in an extended current 
algebra. This current algebra is dual to a doubled Chern-Simons theory in the 3d bulk, 
with one Chern-Simons gauge field for the holomorphic and anti-holomorphic part each (see 
e.g. [29]). Thus the action (2.1) may appear somewhat unfamiliar, and though it has been 
studied in a number of recent works its correct holographic interpretation has been debated 
[30]. 

We claim that the action (2.1) can actually be holographically interpreted as the dual to 
a two-dimensional field theory with a dynamical conserved current. However the theory is 
not conformal: it has a marginal deformation, which is simply j^ja and which must be dealt 
with appropriately. We will discuss the appropriate boundary conditions below: however 
first we will use scalar/ vector duality in 3d to reformulate the bulk dynamics in terms of a 
scalar G rather than the gauge field This is very convenient for dealing with monopoles, 
which are simply point sources for the field B. 



A. Dualizing gauge fields in the bulk and AdS/CFT 

We begin with the bulk gauge field action in Euclidean signature 

Sa = J d'X^^F' . (2.5) 

It is well-known that in three bulk dimensions the dynamics of a [/(I) gauge field is identical 
to that of a scalar via dA ^ ^rfB. To see this at the level of the path integral we note that 
it is implicit in (2.5) that the dynamical variable is the gauge field A, i.e. 

Zbuik = J[DA]exp{-SA) . (2.6) 

Consider replacing the integral over A with one over the field strength F instead. In that 
case the Bianchi identity dF = will no longer be explicitly satisfied, and so we should 
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introduce a Lagrange multiplier field B to enforce it: 

Zbuik = j [VF][VQ] exp (^-Sa + d'X^ e(X)e^'^^a^F,,(X)^ (2.7) 

The normalization of G is arbitrary at this point: we shall explain the rationale for the 
factor of i later. The action for F is quadratic, and we can integrate it out by imposing its 
equation of motion in the action, 
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after which we find simply 



Zbuik = J [ve] exp (^-^ J d^'x^ivef^ , (2.9) 

i.e. the action of a free massless scalar. Note that the dynamical Maxwell equation rf^F = 
has become an identity in terms of 0, while the Bianchi identity dF = is equivalent to the 
dynamical scalar equation of motion di^dO = 0. 

The preceding manipulations are entirely standard. Now we should interpret them from 
the point of view of AdS/CFT. Usually we interpret the boundary value of the gauge field 
cxo) as determining the source for the boundary theory current j^^ i.e 

To write this boundary condition in terms of the bulk scalar we note that a conserved current 
in two dimensions can be related to a scalar operator O as 

f = ie''^dbO . (2.11) 

where the i is natural if we are working in Euclidean signature. We insert this on the 
right-hand side of (2.10) and integrate by parts, after which the source term becomes 

jWS^^O. (2.12) 

Thus the operator O naturally couples to the field strength of a, which in two dimensions 
is a scalar. However in AdS/CFT the field strength of a is precisely the boundary value of 
the bulk field strength F^iy(r cx)), which from (2.8) is: 

'-e-'Fa, = gl^g^^drQ (2.13) 

The right hand side is the bulk canonical momentum 11 conjugate to B with respect to an 
r-foliation. Thus in the scalar representation it is the boundary value of 11 that is the source 
for the field theory operator O. Note that a massless scalar B in AdSs has two possible 
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quantizations with different dimensions of the boundary theory operator A± = 1 ± 1 [31]; 
choosing 11 as the source picks the dimension A_ = 0, which is indeed the dimension of 
O. Note that duahzing a gauge field in 3 + 1 bulk dimensions has a similar action on the 
boundary conditions of the dual gauge field [32, 33]. 

Thus to compute two-point functions of the current from AdS/CFT in the scalar 
representation, we should use the standard holographic prescription to compute boundary 
correlation functions of the operator O dual to B with the quantization appropriate to A_ 
and then use 

(/(yi)/(y2)) = g^^d^A<^(y^Wy^))- (2-14) 

There is thus only one nontrivial function's worth of information in the current-current 
correlator: in two dimensions current conservation strongly constrains the correlator. Note 
also that using the standard AdS/CFT expression for the boundary theory current 

if) = ^VgF^'-ir ^ oo), (2.15) 
9f 

we can express the current in terms of the scalar using (2.8): 

(2.16) 

Thus local conservation of the current is now an identity. Comparing this to (2.11) wee see 
that {O) = B|^^oo which makes it clear we are allowing B to fiuctuate on the boundary. 

In Lorentzian signature we have (j^) = dxQ{r oc). If we consider defining the field 
theory on a circle by identifying x ^ x + Lx^ then the total field theory charge Q = J^^ dx{f) 
is simply the total winding in B at infinity. In particular, the requirement that the field 
theory charge be quantized in units of means that B should be a periodic variable, i.e. 
B - B + ge. 



B. Marginal deformations and boundary conditions 

We have now understood a scalar representation of pure Maxwell theory on AdSa. We 
previously mentioned this theory can be understood as being holographically dual to a two- 
dimensional field theory with a conserved current and a marginal deformation j^ja- This 
deformation is not exactly marginal; depending on its sign it is either marginally relevant 
or irrelevant, resulting in logarithms in various expressions. We now demonstrate these 
facts with a systematic study of the appropriate boundary conditions on bulk fields. For 
convenience we use the scalar representation derived above. 

For this section we will work on pure Euclidean AdSa with metric (2.3). It will be helpful 
to work with the canonical momentum 11 conjugate to B, which we take to be 

n = -^g^'drS . (2.17) 
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We work in Fourier space with 



e(r, X, r) = e(r)e^^^+^^^ / = cj^ + k\ (2.18) 
and the massless scalar equation of motion is then 

2 

dr{r^dre) -^e = 0. (2.19) 
A near-boundary analysis of the solutions to this equation yields: 

e{r ^ oo) ^ A + ^ (^B - A^logr^ . (2.20) 

There are two possible boundary conditions on 0. In one of them, we take the boundary 
value of B(r ^ oc) itself to be the source, i.e. standard quantization. This is a conformally 
invariant boundary condition, as the logarithmic term in (2.20) is subleading at large r. 
It corresponds to the dual operator O having a conformal dimension A+ = 2. Following 
the discussion of the previous section this quantization does not correspond to a dynamical 
current. Indeed from (2.16) we see that fixing the boundary value of the scalar corresponds to 
fixing the field-theory current, and thus this quantization corresponds to having a dynamical 
U{1) gauge field on the boundary, as is discussed in [30, 33]. 

To obtain a dynamical current we thus have to choose the other quantization, corre- 
sponding loosely to fixing the value of the canonical momentum 11 at infinity, which via 
(2.13) is equivalent to specifying the boundary electric field Fxr that the fiuctuating current 
feels. This corresponds to the "conformal" dimension A_ = 0. However the logarithmic 
term in (2.20) means that 11 at infinity actually changes slowly with r, never approaching a 
constant. Thus there is an ambiguity in what we mean by "the boundary value of 11" . From 
(2.20) at infinity we have 

--^n(r) = -/©(r) (2.21) 
dlogr 

Thus any boundary condition cannot be imposed at infinity and must be imposed at some 
specific value of r that defines a cutoff, say r = r\. Consider then the following general 
boundary condition for a normalizable perturbation: 

,^U{rA)=p'e{rA) . (2.22) 

Note that we are fixing a linear combination of 11 and G at infinity, and is thus a double- 
trace coupling, which runs logarithmically [34]. To understand this consider instead imposing 
the boundary condition at a different value of = Xr\ with a different value for the coupling 
Then from the expansion we have 11 (r^) = 11 (rA) — Qp^ log A, and thus the new boundary 
condition can be written 

-n(rA) = /e(rA) . (2.23) 



1 + k' log A 
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This is actually equivalent to the original boundary condition (2.22) if and are related 
by 

^ = --logA. (2.24) 

This is precisely the behavior of a logarithmically running coupling. 

To have a well-defined problem we specify the value of k at some scale; to understand the 
physics at a different scale we can use the evolution law given by (2.24). As usual, physical 
observables should not depend on and Ta individually: rather, dimensional transmutation 
occurs, and all physical observables can be expressed in terms of the RG-invariant scale 

n = r^e^ . (2.25) 

This scale should be viewed as the strong-coupling scale of the problem. Note for > we 
have > ta, and the coupling is marginally irrelevant; if we study only infrared observables 
then its eflFect is benign, only introducing logarithms in various expressions. This is the case 
we will study in this paper. 

If < 0, then < ta and this coupling is growing strong in the infrared, probably 
resulting in new infrared physics. While we do expect this to be interesting, we will not 
study this case further in this work. Note that = separates two theories with very 
different infrared physics and so formally corresponds to a quantum critical point, more 
specifically one of the type called "marginal" in [35] . 

We now interpret the boundary condition (2.22) in terms of the conserved current. Using 
(2.16) and (2.13) we have 

^^'da{j,) = ^^e^'F,,{r^). (2.26) 

This expression may appear familiar: the left-hand side is the divergence of the axial current 
^ahj^ in the field theory, and Fah is the field-theory source. In a 2d conformal field theory 
an expression such as this is true at all scales, with 4- replaced by a precise integer, the 
level of the current algebra k. Indeed this expression would have been the axial anomaly 
in 2d, and should be thought of as a Ward identity that allows a current to be split 
into separately conserved holomorphic and antiholomorphic parts. In our theory, however, 
this expression is only true at a particular r-value. This is dual to the fact that requires 
renormalization, spoiling the current algebra structure. This logarithmic running of is the 
crucial difference between Maxwell electromagnetism in the bulk (which is dual to only a 
single conserved current) and a doubled Chern-Simons theory (which is dual to a Kac-Moody 
current algebra) . We will comment further on the relation between these two theories in the 
conclusion; for now we will simply study the Maxwell theory. 

Note that a boundary action which is consistent the boundary condition (2.22) can be 
constructed and is given in Appendix B. It contains a term which we can directly interpret 
as the current double trace coupling: 

s. = -f I dMyef = f f ff , (2.27) 
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where we have used (2.16). 

For future reference we note that if want to apply a field theory source H{y) for the 
operator O defined in (2.11) then the appropriate boundary condition is (from (2.22)): 

n(rA, y) + -V'e{rA, y) = H{y), (2.28) 

where y only runs over field theory directions. 

Using (2.14) we conclude that the current-current correlation function in momentum 
space is: 

1 e^^e^^PcPd 



{f{p)f{-p)) = 



(2.29) 



9l (n/e)-ip2|_^ 

For illustrative purposes we compute this two-point function for pure AdSs (2.3). The 
solution that is regular in the interior is a Bessel function B(r) = r~^Ki (^); expanding this 
at infinity we find 

{fip)j\-p)) = ^i!!^!^ = 2ne-\ (2.30) 



3f log 



As claimed ta and have combined (modulo factors of 2 and the Euler-Gamma number) 
into the RG-invariant combination r^. The correlator runs logarithmically with momentum: 
a CFT correlator takes exactly the same form except with the logarithm replaced by a 
constant. 

Note that the correlator appears to have a pole at = f^. From (2.25) for an irrelevant 
coupling > this pole lies at ^ r\: physically this pole is at a momentum far beyond 
the cutofl[ and is spurious. In fact even the mathematical analysis leading to it is incorrect, 
as the expansion of the Bessel function used to obtain (2.30) is not valid for p ^ r\. Such 
singularities, which turn out to be ghosts, were discussed carefully in [39] and their removal 
by an appropriate UV completion is discussed in [40]. 

C. Charged BTZ black hole 

We will now write down a black hole solution of (2.1) with an electric fiux emanating from 
the horizon. This is the charged BTZ black hole [24, 25]. Similar states have been studied 
in the context of applications of holography to condensed matter in [41-43]. It defines a 
field theory state at nonzero charge density, which will be the subject of our investigation. 
The metric and gauge field are: 



ds' = f{r)dT' + r'dx' + , A = i( -f^ plog(r/r,)rfr , (2.31) 

}[r) \IQ7iGnJ 

where /(r) = — — [J? j'i) log(r/r/j). The factor of i in A is there because we are working 
in Euclidean. It is useful to think of p as the energy scale in the field theory below which 
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the charge density becomes dynamicaUy important. The real charge density is related to 
this parameter by a large factor (of order 'W^"): 

(2.32) 

We can also write the solution in terms of the dual scalar where non-zero charge density 
imposes a boundary condition on the winding of Q{x + Lx) = B(x) + pLx around the x 
direction. Using (2.8) the required solution is then simply: 

e = px (2.33) 

This fact - that a radial electric field corresponds to a dual scalar that depends on the 
direction parallel to the horizon - is the ultimate origin of the Friedel oscillation structure 
in the boundary field theory correlators. 

The black hole thermodynamics can be worked out. The horizon lies at r = r/^ and the 
Hawking temperature and chemical potential are: 



^1 - = (t7^ ) plog(r./r,) . (2.34) 



In computing the chemical potential we have included a contribution, on top of the usual 
AdS/CFT answer, from the double trace current deformation given in (2.27): 

^^ = A{r^) + ^/^(//) = A(rA) + (2.35) 

As expected the two contributions combine to give the RG invariant log(r^/r/^) in (2.34). 
The zero temperature limit reveals an extremal horizon at r/^ = 2p with non-zero entropy 
density S (x p/Gn and an associated AdS2 x K near horizon geometry. These specific 
features of the zero-T limit will not be important to our results which will continue to hold 
in other more general backgrounds with charged horizons: for example one could consider 
three dimensional versions of the backgrounds discussed in [36-38] . 

More importantly, there are various instabilities associated with the horizon radius be- 
coming too large compared to the UV scale r^. For example at T = the compressibility is 

dp ^ I 1 ^2 36) 

dp gllog{r^/rh) - I' 

which diverges at Th = r^/e and is negative for larger r/^. A negative susceptibility indicates 
a thermodynamic instability. This is not surprising: we should UV complete this theory at 
some energy /radius below the scale and we do not expect black holes to make any sense 
in our current setup beyond this radius. We note that this is also an obstruction to making 
sense of the theory when is an IR scale perhaps associated with some strong dynamics. 
That is when < such that < ta from (2.25). As discussed above, in this paper we will 
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only study the case when is a UV scale, and so we will always have ^ Vh with stable 
thermodynamics. 

Finally, if we consider a probe limit by taking Gn ^ ^ while holding g^p fixed, then p ^ 0. 
As the theory is conformal (modulo logarithms) the dimensionless parameter characterizing 
the temperature is ?, which goes to infinity in this limit. Thus, as is usual in AdS/CFT, 
the probe limit corresponds to a high temperature limit. We will work in this probe limit 
in Section IV. 

III. MONOPOLE SUMS IN ADS 

We now describe the nonperturbative objects of interest in this paper: magnetic 
monopoles. We first review some basic features, and then we build the framework for 
summing over monopoles in AdS and extracting corrections to boundary theory correlation 
functions. 

A. Monopole generalities 

In three Euclidean spacetime dimensions a monopole is a localized instanton that creates 
magnetic fiux. Near the location of the monopole at we have 

e^''Pd,F,,= '^-^5^^\X-Xr^), (3.1) 
V 9 

where is the magnetic charge of a single monopole. Monopoles exist only in compact 
[/(I) gauge theories, i.e. they must satisfy the Dirac quantization condition 

?e?m = 2n7r, n G Z (3.2) 

with the unit of electric charge in the theory. These instantons correspond to tunneling 
events in which the total amount of magnetic fiux in the system 

$= / f (3.3) 

(with S a spacelike slice) is changed by q^- It is clear from (3.1) that the description of 
the field strength in terms of a gauge field F = dA is breaking down at the location of the 
monopole; however the scalar representation introduced in Section II A is ideal for treating 
monopoles, as we now review. 

We begin by considering the partition function Zi^uik in the presence of a single monopole 
at Xj^. Recall the scalar B was introduced as a Lagrange multiplier to enforce the Bianchi 
identity. However we no longer want dF = everywhere; rather we need to incorporate the 
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monopole source as in (3.1). (2.7) is then modified to read 

Zbuik[Xj^ = j [VF] [VQ] exp(^-SA-S,+ '-J d'X^ e{X) {e'''''^d^F,,{X) (3.4) 

-^-^S^'\X - XJ)^ (3.5) 



where Sc is understood to be the action cost associated with the core of the monopole. 
Ehminating F we find 

Zbuik[X^] = J [Ve] exp (- Y y d^'X^iVef -S,- ^g^e(X^)) , (3.6) 

The last term is the coupling of the field to a heavy magnetic source located at a specified 
point Xjri. Note this is the magnetic equivalent of the coupling of the gauge field to a heavy 
electric source iqe A located along a specified worldline C. The dynamical equation of 
motion following from varying this action 

glV'eiX) = tqj-!-^^^^^^ (3.7) 

is equivalent to (3.1). This action will be our starting point for subsequent calculations. 

It is convenient at this point to introduce the bulk-to-bulk scalar propagator G, which 
satisfies 

V^GiX, X') = —5(3) {X - X') . (3.8) 

The field measured at X produced by a monopole at is Qxm{^) = ^-p^G{X^ Xm)- We 
will construct this propagator explicitly on various geometries in future sections. 

The preceding considerations are entirely standard. However we are interested in studying 
monopoles not in vacuum but on the charged black hole background (2.31). We have a 
background electric field and so the field B has a background value Bq = px, i.e. it is 
linearly increasing in the direction parallel to the horizon. The source term in (3.6) now 
means each monopole contributes an extra term to the action. Evaluating the source term 
on the background, we find 

Sberry[Xm] -^<?m©o(^m) -iqmpXm • (3.9) 

This phase acquired by a magnetic charge in an electric field is entirely analogous to the 
usual Aharonov-Bohm phase accrued by an electric charge moving in a magnetic field, which 
is found by evaluating iq^J^A along the trajectory of the charge. In the monopole sum, 
each monopole event should be weighted by such a phase in the functional integral; each 
monopole truly knows where it is along the x direction. We discuss some elementary aspects 
related to this phase in Appendix A. 
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B. On-shell action and Witten diagram 



Now we perform the sum over monopoles. Before starting we note that we assume 
throughout that the field B has a background value 9o = px from the background electric 
field. The calculations in this section will concern perturbations around this background, 
i.e. we split 

e = eo + ^(x) (3.10) 

and deal with the perturbation field 0{X). As Bq solves the background equations of motion 
the action for the perturbation is essentially quadratic, except for a linear contribution from 
the Berry phase term (3.9) at the location of the monopole sources, which we will explicitly 
include. 

Thus our starting point is the action for a configuration of a total of monopoles and 
N anti-monopoles located at points {Xi} (3.6): 



= j d'X^ (y^'^oA +J2(Sc± iqmOiXi) + S,erry[X,]) (3.11) 

Xi 

Here the choice of sign in ± depends on whether we have a monopole or anti-monopole at 
the corresponding point. Finally, to enforce a good variational principle the action requires 
boundary terms that we have not explicitly written down; these will be discussed below. 
The full bulk partition function can be written as a sum over different monopole sectors, i.e. 

Z = z(O'O) + Z(i'O) + + Z(i'i) + . . . (3.12) 

The term Z^^'^^ has no monopoles and is the usual perturbative AdS/CFT answer. Here 
we will compute Z^^'^\ corresponding to a single monopole-anti-monopole pair. This is the 
leading correction; terms with net monopole charge will not contribute to this particular 
observable. We have 

Z(i'i) = j [Dd]K^e-^'^ j d'X^d'X^^giXj ^/ g{X^) exp[-5[X^, X^]] (3.13) 

We place the monopole at A^ and the anti-monopole at A^. To determine the full partition 
function we integrate over their locations, assuming that the correct measure is simply the 
proper volume with respect to the bulk metric, with the scale provided by A, which we 
take to be a hulk UV scale associated with the core of the monopole. Finally, as always in 
AdS/CFT we need to specify boundary conditions on the fields. We allow for a nontrivial 
source H{y)^ and the precise boundary condition at infinity is (2.28). 

We will now evaluate this functional integral by saddle-point. We should thus determine 
the on-shell action on a monopole configuration. The equation of motion for the field 6 is 

9ly'0 = Y,[±^^5^^\X-X,)^ . (3.14) 
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Integrating the action by parts and evaluating it on-shell we find for the total classical action: 



Sonsheu[{Xi}] = J] (^^c ± ^^^^(^i)) + S,,rry[{Xi}] + Sd[e{r^)l (3.15) 

Xi 

where Sd[6] is a boundary term whose form we will describe below. 

Now we would like to compute the current-current correlator in the boundary theory. Via 
(2.14) this is simply related to the correlator of the operator O dual to the field B, whose 
two-point function is 

We are only interested in the contribution from Z^^'^\ Putting the on-shell action into (3.13) 
we find that the correction to the correlator yields 



^2 



SH{y,)SH{y2) 



~2^m (^(-^m) — — Sherry — Sq 



(3.17) 

For notational convenience we will not write down the partition function in the denominator 
Z[0]; as usual its effect is to remove disconnected diagrams. Thus essentially we evaluate 
the field on the locations of the monopoles and integrate their locations throughout the bulk 
spacetime. We now need to determine 0{X). This is the solution of (3.14) subject to the 
appropriate boundary condition and can be written formally as 

e{X) = ^ {GiX, - G{X, X^)) + / £y'K{X; y')H{y') . (3.18) 

Here G{X^ X') is the bulk-to-bulk propagator (3.8) constructed in the previous section. 
K{X]y') is the bulk-to-boundary propagator that satisfies the boundary conditions (2.28). 

Now we put this into (3.17). There are several terms. We first focus on the terms coming 
from the bulk-to-bulk propagators, which are 

^hh = (2G(X^, X^) — G(X^, X^) — G(X^, X^)) (3.19) 

Here the first term is a long-range interaction between the two monopoles. Each of the 
second two terms is clearly singular, corresponding to the infinite self-energy of a point 
charge. One might have thought this self-energy would be included in the single-monopole 
action Sc\ however the dependence of these singular functions on the holographic coordinate 
r contains extra information, which can be thought of as the energy of the interaction of the 
monopole field with the nontrivial geometry and boundary conditions. In the calculation it 
appears as a contribution to the monopole fugacity, which now depends on r. We discuss 
this in the next subsection, and for now we simply call it Sseif{j\ making the full action of 
the monopole 5'^(r) = S'c + Sseif{r). 
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We now turn to the boundary term Sq. It is shown in the Appendix in (B7) that on the 
monopole solution this takes the form 

So = lJ^ d'y H{y) (^iqUK{X^; y) - y)) + gj. J d'y^K{r = ta, y^; y)H{y') 

'(3.20) 

The last term knows nothing of the monopoles; it is essentially the tree-level answer, and if 
we evaluate it on a monopole configuration it results in a disconnected diagram that should 
be dropped. However the other terms must be included. 

Collecting all of these terms we see that the full quantity being exponentiated is 



exp 



Sm(j'm) 'S'm(^m) ipQm{,XYn -^m) 2~^("^™' ■^Tn} 

9f 



iQm j d^y' {K{X^; y') - K{X^; y')) H{y) 



(3.21) 



Note that in the last term half of the contribution comes from So and half from the contribu- 
tion localized on the monopoles in (3.15). Putting this into (3.17) and taking the functional 
derivatives we find the following expression for the boundary theory Green's function: 

G{yi,y2) = -ql^ j d^Xm,fn^/g^^m ^/g^^fn{K {Xm] yi) " K {Xjn] yi)){K [X^] y2) - K{Xrn]y2))x 

2 



exp ( -ipqm{Xm - Xfn) - %G(X^, Xrn) 

9f 



(3.22) 



where we have defined the radially dependent monopole fugacity = A^e"'^'^^^*^. Each 
term in this expression can be represented as a Witt en diagram as in Figure 2; the two 
boundary insertions are connected by bulk-to-boundary propagators to monopoles, and these 
monopoles are then integrated over the bulk with an exponentially suppressed weight that 
takes into account their interaction. There are four terms in the expression above, cor- 
responding to different ways to attach the boundary points to the monopoles, as shown in 
Figure 3. We expect nontrivial correlations only when the two boundary points are attached 
to diflFerent monopoles. We now compute one such diagram (e.g. diagram A), denoting it 
by ^(?/i,2/2)a- The other interesting diagram (diagram B) is obtained by switching yi and 
?/2. 

It is convenient to switch to momentum space p = {(j^E^k) in the boundary coordinates 
2/1^2, after which we find 

G{p\p^)a = -ql, J d^Xrn,m^g^^'rn,rnfp'{^rn)fp^{rfn) exp (^plym'' + wlyfn'' 

q^ \ 

- ^pqm{Xm - Xrn) - -^G{X^, X^) , (3.23) 

9f j 
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FIG. 2: Witten diagram illustrating one of four terms (the one we compute) in (3.22). 



where the /'s are the mode expansions of the bulk-to-boundary propagators, which should 
be understood as connecting a boundary excitation with momentum p to a bulk point (r, y^): 

K{r, y'^-p) = exp {ipaV'') fp{r) (3.24) 

We now perform the bulk integrals over the monopole locations. It is convenient to define 
a relative and a center of mass coordinate in the field theory directions, i.e. 

^ c - ^ (3.25) 

The integral over can now be done explicitly to give a delta function, as there is 
no dependence on the center of mass coordinate. Thus the answer is now G{p^^P^)a = 
{2^f5^^\p^+p^)g{pi), where 



(7^ 

zp«A^ - ipqmK - ^G{rm, rjn] A) 
9f 



Q{P) = -(L J drrn,m^ygm,m^L,fnd^^fp{^rn)f-p{rfn) exp 

(3.26) 

This is as far as we can go in the abstract. To proceed we need an explicit expression for the 
potential energy G(ri, r2] A) of the two monopoles. This depends on the detailed geometry 
(and indeed, in some cases, on the fluctuations of the geometry sourced by the monopoles), 
and so we will treat the high and zero temperature cases separately. 



IV. HIGH TEMPERATURE 

We first consider the case when the temperature is very large compared to energy scale 
associated with the charge density, ? ^ (X). As discussed earlier this is a probe limit: the 
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FIG. 3: Different contributions to final answer (3.22), corresponding to different ways to attach 
boundary points to the bulk monopole. We only compute A; B is related to it by switching yi and 
^2- C and D are not expected to contribute interesting correlations. £^ is a disconnected diagram 
that can be seen in (3.20): it is essentially equivalent to the tree-level answer and will cancel if we 
compute a properly normalized two-point function. 

charge density does not strongly affect the dynamics, and perturbations of the gauge field 
decouple from those of the metric. This greatly simplifies the calculation: we are now dealing 
only with a fluctuating scalar on a fixed background. Thus we take the background to be 
the limit of (2.31), 

i.e. the ordinary uncharged BTZ black hole in planar coordinates. Note that ^^^G ^ and 
so we are still working at a nonzero charge density: it simply does not backreact on the 
metric in the high temperature limit in which we are working. We will often work in Fourier 
space for the perturbations with 0(r, x,r) = ^(r)e^^^'^+^^^, where uje = lnTn with n G Z 
are the Matsubara frequencies. 

A. Single monopole field 

We first solve for the field of a single monopole on this background. We seek a solution 
to the equation (3.8), reproduced here: 

V^GiX, X') = —S^^\X - . (4.2) 
V9 
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where the field 0(X) due to a monopole at point is then 9x^(X) = X^). This 

is just the usual equation for the bulk-bulk Euclidean propagator, and we can use standard 
techniques to solve it. In momentum space the propagator can be expressed in terms of 
solutions to the homogenous wave equation: 

G{r^r;^.,^k) = ^^^^^ (4.3) 

with ^5 the solution that is normalizable at infinity, 9i the solution that is regular at the 
Euclidean horizon, r> and r< respectively the larger and smaller of r and r', and W the 
Wronskian: 

W[a, b] = V^/'-((a,a)5 - idrb)a) (4.4) 
We now construct the solutions 9^ and 9i. The explicit homogenous wave equation is 

dririr' - rl)e\r)) - (-^u^l + e{r) = (4.5) 

As the BTZ metric is locally the same as AdSa, we have enough symmetries to solve the 
wave equation explicitly for all frequency and momentum. The solution that is regular at 
the Euclidean horizon is 

2 \ / 0-2 \ 47rT 



AttT ' AttT ' 27tT 



^^(0 = 2F,(^^,^^;1 + ^;1--^ , (4.6) 



W = 9,{U,- ^ft 



where 2F1 is the hypergeometric function. Now the Wronskian is independent of r and so 
can be evaluated anywhere. A convenient place to evaluate it is at r = ta, where we use the 
fact that 9i) satisfies the boundary condition (2.22) to find 

p^ = ujl + e . (4.7) 

Using (4.6) it is possible to evaluate this Wronskian for all uj and k. However we will 
eventually be interested in the monopole field at large distances, which will be dominated 
by the uje = Matsubara mode and small k. In this limit if we normalize 9i{r ^ 00) = 1 
at infinity, the Wronskian simplifies to 

WiicjE = 0, A; ^ 0) = 9t{rA)e log (^^^ (4.8) 

Note that as claimed r\ and tv have combined into the RG-invariant combination = r\e^ . 
The fact that the Wronskian vanishes as ^ means that there is a pole in the propagator 
and hence a a long range potential between monopoles separated in the x direction. To 
determine the coefficient of the pole we note that in the limit uje = and k ^ both 9i{r) 
and 9i,{r) coincide and become constant, meaning that we have simply 

G(r, iujE = 0,k^0) ^ (4.9) 

Flog(^) 



22 



In position space we then have 

^ r T\r\ 

G{r,r';T,x^oc)=Tj2 / ;^g(r, r^; fc)e-^-+-^- ^ ' (4.10) 

tT^i 27r 21og UBf 



where the last expression correctly captures the leading dependence at large x. We note 
that the problem has effectively become one-dimensional, as in Figure 4. This should not 
be surprising; essentially the field lines cannot spread out along the Euclidean time circle 
as it is compact, and the boundary conditions imposed by the AdS kinematics do not allow 
them to spread out in the AdS radial direction either. Thus they can only extend in the x 
direction, and the linearly growing potential in x is standard for a point charge living in one 
flat dimension. 




FIG. 4: Monopole-antimonopole pair at large spatial separation and finite temperature. Field lines 
cannot spread out in r (due to compactness) or in r (due to AdS kinematics), making the problem 
essentially one dimensional. 



Note also that the sign of the potential depends on whether is an IR scale or a UV 
scale. If it is far larger than the temperature - as we have been assuming - then the overall 
sign is positive and two monopoles of opposite charge will feel an attractive force, which is 
physically reasonable. 



B. Priedel oscillations at high temperature 



Recall that at the end of the previous section we had established the form of the instanton 
correction to the current-current correlator to be (3.26): 

,2 



iPaA"" - ipqm^^ - ^G{rm, Trn] A) 



G{p) = -Qm J drrn,m^gmML,md^^fp{^rn)f-p{rfn) exp 

(4.11) 

With the above expression for the monopole-monopole correlation we are now in a position 
to finish the computation of this correction. We expect the answer to be dominated by the 
contribution from monopoles that are widely separated in the x direction. In this regime we 
can directly use (4.10), G(r,r'; A^, ^ oc) = T^. \A^\. 

At this point we also set the external frequency cj^ = 0, keeping only the spatial momen- 
tum k. In this limit at large A^^ nothing depends on A^-, and so we can explicitly do the 
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integral, picking up a factor of the inverse temperature. The remaining integral over A^^ is: 
I{k) = - / c^A.exp (^A:A, - ipq^/\^ - 2nTK^\/\^\) = ^ (4.12) 

where the definition of the parameter is for later convenience. This expression should 
be understood as correctly capturing the contribution to the integral from large A^;. The 
integral is a Lorentzian and we find 

I(k) = 4:7rK^ ^ ^ (4.13) 

This is the interesting momentum-dependence of the final answer. The peak of the 
Lorentzian is at a A; value that is shifted by pg^; this is precisely the desired Friedel os- 
cillation, and the wavevector of the oscillation could clearly have been predicted from the 
Berry phase term alone, with essentially no computation. However the detailed form of 
the finite-temperature Lorentzian followed from detailed considerations involving the bulk 
geometry. We leave further discussion of the result to the conclusion section. 

We now note that there still remains a momentum-dependent radial integral to do. The 
full answer takes the form 

G{k) = -qlK^R{kfl{k) (4.14) 
where the radial integral has been absorbed into a function, 

rfrV^e-^-MA(r) (4.15) 

and we have used the fact that R{k) = R{—k). This expression determines the eflFective 
contribution of a monopole at radius r to the answer. In the field theory each monopole 
corresponds to some sort of instanton event, with its radial position r loosely corresponding 
to the size of the corresponding instanton. Thus this expression can be interpreted from the 
field theory point of view as some sort of monopole "density of states" . 

//e(r) is the mode expansion of the bulk-to-boundary propagator, and so is simply the 
(appropriately normalized) solution to the wave equation that is regular at the horizon (4.6). 
We are interested in this expression at very large momentum, k ^ q^p ^ N'^- Note that 
this very large Euclidean momentum means that we are very off-shell, and we expect the 
mode function to fall off quickly as we penetrate to the bulk of the spacetime. While this 
can be confirmed from examining (4.6), it is simpler to realize that at high momentum we 
are probing the high-energy CFT2 structure of the theory, and thus (at least at r Th) the 
mode function should be the appropriate solution to the wave function on pure AdSs in the 
Poincare patch, i.e. the Bessel function 

K (-) 

fk^ooir > Th) '-^^ . = 2r,e-^ (4.16) 

rklog 
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Here the normalization has been fixed by requiring that at infinity fp{r) corresponds via 
(2.28) to a source of unit magnitude for aU momentum: H{k) = 1. Note that in order to 
avoid problems with the UV pole at A; = we must demand that k ^ q^p <C r^. This is 
easy to achieve with a moderately small ^ l/log(A^). 

Note now that at intermediate radii r/^ <C r <C we have from the expansion of the Bessel 
function //c(r) ^ exp (— ^) , and thus the mode function is indeed exponentially suppressed in 
the bulk. This means that (modulo contributions from Sgeifi^)) the bulk of the contribution 
comes from monopoles very close to the AdS boundary; the corresponding field-theoretical 
instantons are almost as small as can be, as the large ones do not couple strongly to the 
hard probe that is a high-momentum correlation function 

There is also some radial dependence in the function ^^^(r). As explained in detail in 
Appendix C, this corresponds to a radially dependent monopole fugacity. While it can in 
principle be computed in terms of a sum over the mode functions, this is numerically very 
time-consuming for the BTZ black hole. We instead compute the fugacity for pure AdSa, 
as the dominant contribution to the monopoles comes from radial positions r ^ and we 
expect this to be a good approximation to the answer. It is 

2 - 

Sm{r) = S^ore " ^ log log f ^) , (4.17) 



This result is discussed further in the Appendix. It wants the monopoles to be towards the 
interior of the geometry, but clearly depends very weakly on r and arises purely from the 
marginally running boundary conditions. 

Thus the gentle dependence on r of the fugacity does not change the fact that the integral 
in (4.15) is UV-divergent. There is a cutoff-dependent piece in the answer: 

R^k) ^ (^!^)%-5'"('-A)log-^ (A^ (4.18) 

We may thus assemble all of the pieces to obtain the full form of the Friedel oscillation at 
finite temperature: 

{Sp{k)Sp{-k)) = e-''-(^- \ rj^^T^^td. ^ + ^ (4-^9) 

{2tiTK^ ) +{k- pqm) 

where we have used (2.14) to relate the answer to a density-density correlation and the overall 
constant A ^ Ang^'^A^K^^ log~^ (jf^^ • We have also added in the contribution from the 
other contributing diagram (i.e. diagram B) in Figure 3; this results in an answer that is 
symmetric in k by creating the corresponding singularity at A; = —pqm- In determining the 
prefactor we have set k pqm', as this expression should be understood only as capturing 
the leading momentum dependence in the vicinity of the pq^- 
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V. ZERO TEMPERATURE 



We turn now to the structure of the Friedel osciUations at zero temperature. These 
osciUations will occur at the same wavevector - pg^ - as the large temperature analysis of 
the previous section. However the form of the singularity should be different. In particular 
in momentum space the Lorentzian (4.19), whose width vanishes as T ^ 0, should become 
a genuine singularity whose form we would like to determine. 



A. Sound mode 

While our final goal is to compute the contribution of monopoles to the density-density 
correlator at zero temperature in the charged BTZ black hole, we start by computing the 
answer in the absence of the monopoles. We go through this in some detail, both because 
it is interesting in its own right and we will use the results heavily later in determining 
the form of the long-range interaction between monopoles. Interestingly we will find a zero 
temperature sound mode similar to the sound mode found in the higher dimensional charged 
black hole [11]. This sound mode will turn out to control the form of the Friedel oscillations. 

Since the charged black hole mixes the scalar field (or the gauge field in the dual repre- 
sentation) with the metric the dynamical problem becomes much more difficult. This sort 
of problem has been addressed in a similar context by [11, 12] and we follow their methods. 
The key is to find master variables that diagonalize this mixing. We will call these modes 
"gauge invariant" for reasons which will be clear. 



1. The gauge invariant variables 

We start by perturbing the charged BTZ background (2.31) using the dual scalar variable 
B. We pick a gauge for the metric fiuctuations suited for holographic computations gra = 
and grr = f{r)-^- 



Note we have rescaled B = (IGnGNgp) B for future convenience. The same rescaling of 
p ^ p was introduced in Section II C. The linearization of the Einstein equations and the 
scalar equation for B defines the dynamical problem we would like to solve. The equations 
will be left in Appendix D since they are not very enlightening. 

It is useful to note that there are three unfixed diffeomorphisms {r ^ r + y/JC -> ^ ~^ 



1/2 



26 



X + ^ r ^ T + under which these fields transform: 

K ^ K + 2drC + -yjC , hi ^ hi + 2d,C + , (5.3) 

where the (C, C O's satisfy the foUowing : 

drC = , drC = -^^-'^^ ' = "-^^rC • (5-4) 

There are 4 independent fields /i^, /i'!^., /i^) each with second order equations of motion 
in the radial direction. Together with the 3 first order constraints from Einsteins equations 
and 3 unfixed diffeomorphisms we should be left with two independent first order degrees 
of freedom (or one second order). This is of course expected since in three dimensions 
gravity does not contribute any propagating modes. In order to cleanly solve this problem 
we should write down two degrees of freedom which are gauge invariant under the left over 
diffeomorphisms (5.3). We choose these as: 

a' = D^e - pfd^hl + ^d^hl + 2pr^drhl , = 2rf9' - 2dl9 + pd^hl , (5.5) 

where Di = — {f'rd^ + 2r^5^). From now on we will work in Fourier space with our con- 
ventions given in (2.18). The dynamical equation for the a's may then be derived from the 
equation of motion (D6-D9) and constraints (D3-D5) by brute force. See Appendix D for a 
general discussion of this. They must close on themselves because we know they are gauge 
invariant and we know there are only two independent such variables: 



Start by imposing a standard regularity condition on the a's in the interior at r = r/^. 
For example at T = this solution behaves as follows close to the horizon: 

4^!^ .)(r-n)'^e-^. (5.7) 

+ 2\uje\p J 

Solving (5.6) at large r we find the following behavior: 

= ar^ + ... , = alogr + /3 + . . . , (5.8) 

where the ratio (3/ a will be fixed by the regularity condition. We impose the following 
boundary conditions on the metric: 

hl^O + 0{r-^) , /i; ^ + 0{r-^) , /i", ^ + 0{r-^) , (5.9) 
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since for now we will not be interest in sourcing the stress tensor in the field theory. These 
boundary conditions are enough to fix the diffeomorphisms given in (5.3) and thus we can in 
principle solve for all the dynamical fields subject to the constraints from Einsteins equations. 
So we now have enough information to extract the boundary behavior of 6 using (5.5) and 
(5.8-5.9) we find, 



^ ^ V V ^ ^(^Vr - (/3H - 1/2 - logr) + . . . j . (5.10) 

This can then be used to extract the current current function (after subtracting the back- 
ground density) using (2.29). 

^' glk^-p^{Pla + \og{r.)) ^^'''^ 

Note that it is trivial to analytically continue this Euclidean answer to real times uje ^ —iu. 
See [44] for a discussion of this. In the higher dimensional case the same problem has been 
solved [11] and an examination of the current-current correlator studied in [11, 13, 45] reveal 
no 2A;i^^-like singularities. This should continue to hold in the lower dimensional case at 
hand, despite certain details being different. We have done some unenlightening numerical 
calculations of (5.11) which seem to confirm this suspicion. 



2. Hydrodynamic limit 

We can miraculously solve (5.6) in the small p limit. This is analogous to the usual finite 
temperature hydrodynamic limit, where instead of requiring the momenta to be small in 
units of the temperature we require they are small in units of p. This allows us to extract 
the zero temperature sound mode. For uj\^k^ <C we have, 

I _2 I = I 2(c.|r^+fcV) p^fc^ I I ^2 I ' l^"!^) 

<Ti = cyf a2 = ^{rf-2f) . (5.13) 

We found this solution with the same methods used to construct Goldstone modes of holo- 
graphically ordered phases [46] . This solution is regular in the interior and so we can use it 
to read off: /3/a = — log{rh) + ^^/p^ — 1/2. The density correlator in p space (5.11) in the 
hydro limit is: 

1 Zk'^ 

WpW-rt> = ^ ^^^„,,,^^ , (5.14) 



which has a solution 



where, 



^ v^^ = l-Z. (5.15) 



\og{r^/rh) + 1/2 - r^/p2 
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Note that = iSp. Clearly there is a sound mode - in real time there is a mode which 
disperses as = ±Vsk. To make sure this is a genuine sound mode we have to show that 
the width from the self energy S is small in the hydrodynamic limit. This is computed in 
Appendix E where we find the width is highly suppressed T>{(jUe ^ ^ - 
At zero temperature the speed of sound satisfies: 

where we have used the thermodynamics of the charged BTZ black hole discussed around 
(2.36). As discussed in Section II C we are only interested in the regime where r/^ <C 
where it is clear that < '^^ < 1. 

An analogous sound mode was found in a probe brane setup [47] where it was suggested 
to be the holographic realization of the zero sound mode of an underlying Fermi surface. 
Thus giving evidence to the fermionic nature of that phase. It seems reasonable to make a 
similar identification in the charged black hole setup studied here and in higher dimensions 
[11]. A more fine grained comparison of this zero sound mode to a Landau Fermi Liquid 
(LFL) was undertaken both for a probe brane [48] and the higher dimensional charged black 
holes [45] with the former behaving more like a LFL. 

The mode we found will be crucial in revealing the form of the singularity. This is 
similar to what one would expect for a Luttinger liquid, where the sound mode is described 
by a massless scalar and the singularities are created by operators such as e^^"^. For 
further discussion of the relation to a Luttinger liquid see the conclusions. All in all, we feel 
this serves as further evidence for the zero-sound interpretation of these modes. 



B. Priedel oscillation computation 



We are now in a position to compute the monopole anti-monopole potential for monopoles 
situated at large r ^ Th- We concentrate on the monopoles living at large r simply because 
for smaller r their contribution to the Friedel oscillation computation is highly suppressed 
by the bulk to boundary propagator.^ In Section IV B we have already argued this to be the 
case at high temperature, and it remains the case here since the bulk to boundary propagator 
at large momentum is insensitive to the density or temperature and will continue to be given 
by (4.16). 

At large r ^ the form of the bulk fields in momentum space is still (5.10) and (5.9) 
(actually here we are also taking r ^ k^uoE)- In particular the effects of mixing between 
9 and the gravitational fiuctuations are suppressed. Thus we can make use of the formula 



^ Rather the appropriate probes of the monopoles which necessarily have large momentum cannot penetrate 
to smaller r. 
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(4.3) for the bulk-bulk propagator: 

G{r,r ;tu!E,k) = 7. — - — , (5.17) 

where the regular solution 6i{r^) is proportional to that given in (5.10) and the boundary 
normalizable solution is: 

^,(r>) (X 1+ (-l + 21og(r./r>)) . (5.18) 

For our purposes here we have kept only the dominant terms at large r^r' ^ k^ujE^Th. The 
wronskian in (5.17) is then simply the AdS^ wronskian for the massless scalar wave equation: 
Vl^[a,6] r^{{dra)b — {drh)a). Plugging into (5.17) and keeping the lowest order terms we 
derive the unsurprising expression: 

1 Z 

G{r, r'] iujE, k) = — — — — — — -— - = — — — -. (5.19) 

k^ - p\/3/a + log(r^)) ouE,k<t:rh uj% + v^k^ 

This is essentially the same (up to normalization) as the boundary greens function for 6 
since the appropriate bulk wave functions are approximately constant at large r. For the 
interaction energy at large separation we only need the small k^ujE limit. Taking the fourier 
transform, and cutting off the k^ujE integrals at small v^k^ + cc;|; > l/Ljj^ we find at large 
separation, 

G(,,,V,,) = j|-log(^l±^)+.... (5.20) 

Note that the IR cutoff Ljr will drop out of the final result due to a similar IR divergence 
that appears in the momentum integrals that define the monopole self energy, as is discussed 
further in Appendix C. 

We can now evaluate (3.26) and integrate over the separation. This gives us the frequency 
and momentum dependence close to k ^ QmP-, i-^. the zero temperature generalization of 
(4.13): 

I{UJE, k) ^ {vl{k - q^pf + c^l)^ , K = . (5.21) 

This demonstrates the form of the zero temperature singularity that we were seeking. We 
have omitted an 0(1) correction to the exponent, as K is very large, scaling like an inverse 
power of gjp. The key difference from the finite temperature case is that now the gapless 
sound mode creating the long-range correlation behaves like a massless scalar in the two 
field-theory directions, resulting in this power-law structure. The radial dependence of the 
fugacity at large r will essentially be the same as for the high temperature answer so the 
radial integrals and determination of the form factor R{k) are precisely as discussed in 
Section IV B. 
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C. Monopoles interacting with 3d gravity 

We have obtained an answer. However the alert reader wiU have noticed that we did 
not carefuUy deal with the interactions of the monopoles with gravity; instead we were 
fortunate that the dominant part of the answer came from a region where these interactions 
are suppressed. In other words we never wrote down the generalization of the monopole 
partition function (3.13) which allows for fluctuations of the metric. 

This is for a good reason: the fully backreacted problem is both technically and concep- 
tually diflicult, forcing us to confront various issues related to monopoles and background 
electric flelds interacting with 3d gravity. For this particular problem it turns out that these 
issues do not aflFect the answer. The resolution was argued for above: the monopoles want 
to live at large r where the interactions with the metric that caused these issues in the flrst 
place are suppressed. ^ Despite this fact, in this section we still go through the argument 
in some detail, both because it helps us understand some aspects of the above analysis and 
because it might be useful in future investigations. We dedicate this section to understand- 
ing various issues of monopoles interacting with 3d gravity, ending with the generalization 
of (3.13) and the computation which leads to (5.21). 



1. Monopole size 

The flrst issue we must address is the back reaction of a monopole on gravity. To do this 
we start in 3rf flat space with no electric fleld. It is natural to rescale B = {IQuGnqD'^^^Q 
as was done in (5.2). A monopole sitting at the origin has the solution: 

d^e = iq^5\X) ^ e^^, (5.22) 



where we introduced the parameter deflned as, 

\ 9f J 

and which has the units of length. The back reaction on gravity is then roughly: 



(5.23) 



d'Sgr^T^r^idefr.^^ ^ ^9^^^ (5.24) 

So this is large when |X| ^ g^. That is, close to the monopole the perturbative flat 
space solution breaks down, the magnetic fleld backreacts on gravity, and we must flnd a 
new solution. Indeed one is available: the wormhole solution of [49] has long-range flelds 



^ Of course the mixing of with gravity is still important since the long range fields probe the full geometry, 
not just the geometry at large r. 
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corresponding to those of a magnetic monopole in flat space, but is actually a Euclidean 
wormhole threaded by magnetic flux and connecting two asymptotically flat regions. In 
gravity there is thus a natural core size for these monopoles of order - this is the size 
of the throat of the above wormhole. In fact there are several ways to complete the linear 
solution with a non-linear solution, such as for example the original SU (2) 't Hooft-Polyakov 
monopole where the size of the core of the monopole will be determined by the ly-boson 
mass scale (as long as this length scale is bigger than q^.) 

We do not wish to be speciflc about the exact form of this bulk UV completion. This is 
possible as long as the core size or q^ is much smaller than the AdS radius L. We can 
then argue that the "linearized" solution is correct outside the core and imagine matching 
onto the non-linear core using the flat space answer. Then the only effect of the bulk UV 
completion is on the action associated with the core Sc introduced in (3.5). Either way we 
will always have the constraint that <C L. Recall that we flxed units by setting L = 1: 
we briefly restore it in this section only for illustrative purposes. 

The linearized monopole solution in the charged BTZ black hole is then constructed 
perturbatively in qm/L^ 

9 = 90 + {qjL)g, + 0{qjLf , (5.25) 
e = Go + {qjL)ei + 0{qjLf , (5.26) 

where go and Bq = px are the charged black hole solutions without monopoles. This 
expansion will break down close to the monopole, but in a controlled way. It is important 
that when we write this perturbative expansion we are working with flelds that vary on fleld 
theory scales all of order {(r/L^), cj^, fc, . . .} ~ p, where we introduced this rescaled density 
in (2.32) and we reproduce it here: 

p={167rGr,9l)'/^p. (5-27) 
The action for the solution including the monopoles will then have an expansion of the form: 

Se = (So + iqm/L)Si + {qjLfS2 + ...). (5.28) 

Note that the 3 leading terms in (5.28) are already evident from the high temperature answer 
(3.15). The flrst term is simply the free energy of the charged black hole and cancels from 
the flnal answer so we ignore it. The second term is the berry phase term and takes this 
form when written in terms of p, that is Si = ip{xm — ^m)- The last term is the monopole 
interaction energy and fugacities. Actually the only difference between the high temperature 
calculation and the zero temperature calculation is the fact that at flrst order the metric 
is non-zero. At second order the metric will be non-zero in both cases. So we still secretly 
required that <C L for the high temperature calculation. 

Writing the action (5.28) in this way makes it clear we are doing a Euclidean saddle 
computation of some path integral which we will specify shortly. Quantum corrections will 
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be down by factors of Gn- 



2. Monopole position 

Naively following the steps in Section III we should be performing a path integral over 
both g and defined by integrating the monopole partition function (3.13) over metrics 
weighted by the usual Einstein-Hilbert action: 

Z(i'i) = j [Dg][De] j d'X^^giXj J d'X^^g{X^) exp {-S[X^, X^]) , (5.29) 

where the action (written in terms of the rescaled variables defined above) is now: 



s[x^, x^] = ^^-^ (^j (fx^ (^-i? - 2 + ^(ve)^) + ^g^(e(X„) - e(x^))) 

Varying this with respect to 5^,0, the equations of motion for the saddle are: 



(5.30) 



R,. - \g,u {R + 2) = \t% , v^v^e = ±iqm5^^\X - Xi), (5.31) 

i=m,fn 

where we have also defined the appropriately rescaled scalar stress tensor: 

= (v.ev.e - ^^..(ve)^) . (5.32) 

We now come to a distressing point. These equations seemingly lead to an inconsistency 
since the stress tensor is not conserved: 

v%t = v2ea,e= ^ ±tq^dMXi). (5.33) 

i=m,fn 

The right hand side of this expression is not obviously zero. This contradicts Einstein's 
equations: gravity must always be coupled to a conserved stress tensor. Upon further 
refiection, this violation of momentum conservation is expected; it is intimately related to 
the Berry phase and can be understood in terms of the momentum carried by the monopole 
event, as is discussed further in Appendix A. Nevertheless, this is a serious problem since 
we cannot proceed if we do not have a consistent set of equations to solve. 

The resolution to this puzzle is extremely simple. We should have also extremized the 
action (5.30) with respect to the position of the monopoles. Said diflFerently we should do the 



^ Note that various fields are (e.g 0i and gi) are purely imaginary on this saddle. We should not be 
distressed by this. In our Euclidean calculation the classical field configuration we are evaluating is 
related to a quantum tunneling amplitude; the interplay between the real and imaginary parts of the 
action of this configuration carries information related to the imaginary Berry phase and so presumably 
is related to quantum mechanical interference between different tunneling paths. 
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integrals over monopole positions in (5.29) within the saddle approximation - thus letting 
the monopoles roam free in the bulk and telling us where they want to sit. Varying with 
respect to the Xi we find, 

a^e(X,)=0, i = m,m. (5.34) 

This solves the contradiction presented in (5.33). The situation is the same as if we had 
coupled a massive particle to gravity. The world line of the particle must be a geodesic: 
otherwise the stress tensor associated with the world line is not conserved. 

Note that the situation is diflFerent from the case without backreaction. The position 
of of a given monopole is no longer a free parameter which we integrate over in the end. 
Rather it depends on the fields B and which in turn depend on the position of the other 
monopoles. We can trace the reason for this to the fact that diffeomorphisms C,^ move 
around the position of the monopole insertions B(X^) B(X^) + C^5^B(X^) such that 
the gravity path integral cannot be simply disentangled from the integrals over the monopole 
collective coordinates. The infinite dimensional field space B has been mixed up with the 
6 dimensional monopole coordinate space. There is no generic way to choose a 6 

dimensional slice through the extended field space. At high temperatures there was a well 
defined way to do this split, and we left the monopole integrals to the end. Here we are forced 
to do the integrals at the same time as the gravitational path integral. These complications 
appear to be related to the usual issues of not having local observables in quantum gravity. 

At this stage we note that it seems hard to solve (5.34) in the perturbative qm/ L expan- 
sion we outlined in (5.25-5.26) since ^^^B is non-zero at zeroth order due to the background 
electric field. Thus it seems that in order to satisfy (5.34) one must work outside the per- 
turbative expansion - eflFectively canceling the zeroth order term with the first order (q^/L) 
interaction term. It might be possible to do this if we move the monopole anti-monopole 
pair very close together where the interaction potential is large: this is physically reasonable, 
essentially the Berry phase effectively binds monopole-antimonopole pairs together and pre- 
vents them from proliferating. However this is not the physics we are looking for; rather 
we seek a contribution to a holographic two-point function, which should depend on physics 
determined by monopoles separated by a large distance. 



3. High momentum bulk to boundary propagator 

There is one further ingredient in our holographic computation: the bulk-to-boundary 
propagator, which as we will see should have been included in the saddle. In the regime of 
interest the momentum fiowing in the propagator is very large: 

p-^2kF = -TTT^^qmP . (5.35) 

Thus this propagator will actually back-react on gravity in a well defined way which we 
explore now. Indeed it appears at exactly the same order as the Berry phase term and thus 
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can resolve the issue of satisfying (5.34). Following [50] we write a "massless worldline" path 
integral representation for the bulk to boundary propagator from a point y on the boundary 
to a bulk point as: 

Here the integral is over all paths X{s) from the boundary to the specified point in the bulk. 
5 is a parameter along the worldline, and the boundary conditions at 5 = fix the path to 
the boundary at 

f(0) = rA, X^{0) = y% (5.37) 

where ta is the UV cutoflF. The boundary conditions at 5 = / fix the path to a point in 
the bulk X^{1) = X!^ = (rm^y^) = (^m^'^m^^m)- We have essentially fixed a gauge for the 
einbein leaving an integral over /. 

The utility of K^^"^ is that it allows us to define the bulk to boundary propagator for a 
general metric g (asymptotically AdS). This would be formally difficult in terms of solving 
the PDE wave equation on some general metric g. Also since the momentum in the world 
line is necessarily large we expect the "geodesic" approximation to the path integral to be 
a good one. 

Because it is not clear what we might mean by a "massless geodesic in Euclidean space" 
we consider a simple example which demonstrates that in general we must consider complex 
geodesies. These can be understood as complex saddles of (5.36), a prescription for which 
can be given in terms of the usual WKB approximation to the wave equation [51]. 

We take the metric to be the charged BTZ metric. The saddle equations for the path 
X{s) follow from the conserved charges associated with translation invariance P and the 
massless condition: 

Note that the equation of motion for / imposes the massless geodesic constraint''. 

Now consider taking the Fourier transform of K with respect to the boundary coordinate: 

UXm) = I (fy e-^P-^y'-y-^K{X^, y) . (5.39) 

Again performing this integral by extremizing with respect to we see that this fixes the 
charges in (5.38) to be purely imaginary E = iuje and P = ik. The appropriate geodesic 
thus follows an imaginary x{s)^f{s) trajectory. On this complex geodesic (5.39) computes 



^ The variation of the single number / imposes the constraint at one point, but it is preserved by the 
dynamical equations of motion and so must now be satisfied everywhere. 
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the WKB approximation to the massless wave equation [51]. For example if we set uje 
the path (5.38) written in terms of x{f) is simply given by 

dx i sgn(/c) 



(5.40) 



where the sign is determined by the requirement that (dr/ds) < 0. Then (5.39) becomes, 

/,(X„).e.p(-,|£-^). (5.41) 



This is the WKB approximation to the massless wave equation in the charged BTZ black 
hole at large k. For r large it reduces to the result we used previously in (4.16) for the 
bulk to boundary propagator, where it is clear we should take the decaying WKB solution 
for regularity in the interior. The complex geodesic is telling us that the bulk to boundary 
propagator is highly suppressed at high momentum as we move into the bulk. 

Having built some confidence in K^^"^ we now include them in our path integral (5.29). 
This path integral will then directly computes the density density correlation function, 

1 



{Sp{yi)Sp{y2)) 



[Vg][Ve] I d^X^^/g{X^)K['\x^,yr) j d'X^^g{X^)Ki'\x^,y,) 
xexp(-5[X^,X^]) + (mf^m), (5.42) 



where we will not keep track of the overall normalization since we will only compute the 
path integral in a saddle approximation anyhow. The notation may get confusing, so we 
urge the reader to refer back to Fig. 2. We use / = 1, 2 to label the two paths Xj and the 
two boundary points yj where the currents are inserted . We use i = m,m to label the 
monopole/anti-monopole positions Xi. For the particular saddle we consider here they are 
correlated as m ^ 1, m ^ 2. So for example below we always sum over these indices: 



<S = {(^,/)} = {(m,l),(m,2)}. 



(5.43) 



The diagram Fig. 3B has the opposite correlation. The fuU set of equations that foUow from 
varying everything {g^,^, Q,Xi,li,Xi} are 



R 



2 ''^ 



7=1,2 



/ 

^0 



dsi- 



dXi 
dsj 



dXi 
dsj 



^v^e= J2 ±iqmSlj\x - Xi) , 



fxi 



dXf dX^ 
dsj dsj 







d^X,) = ±2ig,,{X,) 



dX^ 



dsi 



dsi dsi 

(and Xiili) = X,) 



(5.44) 
(5.45) 

(5.46) 
(5.47) 
(5.48) 



si=li 
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In (5.45) and (5.48) the + refers to the couphng the monopole at and the — to the 
antimonopole at X^. It was also natural to rescale (s/,//) = {IGttGn /qm){si/h) similar to 
the rescaling of B ^ B. The two paths are fixed to the boundary at f/(0) = ta, yj{0) = yj. 
It can be checked that the full stress tensor coupling to gravity is conserved when all the 
above equations are satisfied. 

The equation (5.48) is the endpoint equation of motion for the position of the path, or 
really the equation of motion resulting from varying the position of the monopoles. This is 
the generalization of (5.34) that we had trouble satisfying before. The missing momentum 
at the location of the monopole is now being supplied by the bulk-to-boundary propagator. 
One can heuristically think of the bulk to boundary paths as strings attached to monopoles 
from the boundary that allow us to pull the monopoles apart. As we shall see, this can 
be achieved only if we probe the monopoles at the correct (i.e. Friedel "2/ci7" wavevector) 
momentum. 

Finally some care must be taken to define the delta function 5^ in (5.44) since the 
path Xj{sj) will generally follow a complex trajectory. The best way to proceed is to solve 
Einstein's equation on complex coordinates which pass through these Xj trajectories. It 
seems possible to understand this in the perturbative expansion of Qm/L outlined above 
and explored further below. Some inspiration from the real time formalism of thermal field 
theory (i.e. the Schwinger-Keldysh contour) will probably be necessary. 



4. Some aspects of the final saddle 



We can solve (5.44-5.48) using the perturbative framework outlined in Section VC 1. In 
addition to the expansions of the field B and g given in (5.25-5.26) we should also write: 

Xiisi) = X'^isi) + 0{qjL) , (5.49) 
= /7 + 0{qjL) , (5.50) 
Xi = X^ + 0{qjL) , (5.51) 

where we will only need the zeroth order term and thus we drop the superscript from now 
on. Then the procedure for constructing the saddle is simple. We start with the zeroth order 
solution, which for B is simply the charged BTZ black hole. The zeroth order geodesies 
are then the same as discussed around (5.38), where the rescaled conserved charges Ej^Pj 
are fixed by the monopole equation of motion (5.48), 

A . / X dfr ^ ^dxj ^ dxj 1 , , 

dsi dsi dri fiy/f{fi) 

Note that the momentum in the geodesic is fixed by the monopole equation of motion to be 
the Friedel wavevector. It will be often convenient to parameterize the worldline in terms 
of f/ instead of //, it should be clear which one we are using. The bulk positions of the 
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monopole can now be found by solving for the the end point of the above paths: 

dr' 

Tiivi) = Ti = Ti , xi{ri) = Xi = xiTi 7== {i,I)eS, (5.53) 

where (r/, Xj) = are the positions of the boundary insertions. The radial coordinate of the 
monopoles f/(//) = is fixed by the final monopole equation of motion, the r component 
of (5.48). 

jj= = ±idrQ{ri, Ti, Xi) . (5.54) 

It should be clear that by symmetry the two are the same = r^. We will come back to 
this equation shortly. Finally to find B and g to first order we plug the zeroth order geodesic 
back into Einstein's equations and the dynamical equations for G and linearize. This will 
then provide sources for the equations that we discussed in Appendix D . Using the paths 
described above the sources simplify to the following: 

(r,^-19,AR + 2)-It^?)' = -f^e{r-rjJ^niniS{T-Tj)Sc{x-Mr)),{5.55) 
(^V^e) = iqmS{r-r^) ^ ± S{t - Ti)Sc{x - Xi) , (5.56) 



i=m,m 



where we have performed the integral over s by expressing it as an integral over f and 
= {dx =F i^dr) is a null vector pointed along the complex geodesic. At this stage if 
we were to attempt to solve these equations we would need to confront the meaning of Sa- 
lt can be understood in terms of a contour in complex x space which passes through both 
monopole paths xj{r) for r > r^. For r < rm the contour can lie on the real axis. Then 
these equations should be solved on this complex contour. 

Fortunately for us, and perhaps unfortunately for the formalism we developed, much 
of the above discussion is not important for the final result. The reason can be traced 
to (5.54), where the right hand side is actually zero to the order in the Qm/L expansion 
we are working at. So it does not seem possible to solve for r^. Rather one can argue 
that = (X) is the solution - the monopoles live at the boundary. In order to probe the 
monopoles - i.e. in order to couple them consistently to gravity - we had to include the 
bulk to boundary propagators carrying a large amount of momentum. Since they carry large 
momentum the effective tension in the paths of the bulk quanta is large and this pulls the 
paths towards the boundary. This is because these propagators are exponentially suppressed 
in Gn upon moving in to the bulk and so the saddle approximation to (5.42) will always want 
to concentrate the monopoles at the boundary. These are exactly the same considerations 
discussed when performing radial integrals in Section IV above. 

Let us now evaluate the answer for the current-current correlator in this framework. 
On-shell in (5.42), we see that the two monopoles in the bulk are actually located as close 
to the radial cutoflF as they can be: furthermore, from (5.53) their positions in the field 
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theory directions must coincide with those of the boundary theory operator insertions, i.e. 
Vm ^ yi^Vrn ^ 1/2- Thus wc cau siuiply evaluate the on-sheU action with the monopoles held 
at this location. 

The calculation of the on-shell action itself was performed in the previous section. We 
can directly compute the correlator in position space, and we find (adding both diagram A 
and diagram B in Figure 3 together) 

{5p{yi)5p{y2)) cos {pqm{xi - X2)) [{xi - X2Y + v^in - rs)^] , (5.57) 

where the prefactor cannot be reliably obtained in the saddle-point framework . This is 
equivalent to (5.21). Note that it may appear we have performed fewer Fourier transforms; 
however as the momentum-space expression is valid only in the neighbourhood of k = pq^ 
in reality the information stored in these expressions is the same. 

We can go a little further than this by allowing the first order correction to B to feature in 
(5.54). By balancing these two terms that appear at different orders in the qm/ L expansion 
we will eventually find a true saddle for the monopole radial coordinates ^ p/ qm- This 

is necessarily in the AdS^ region of the geometry. 

To see this we again argue that at large r mixing with gravity is suppressed (by factors 
of order p/r.) Additionally the worldline sources in (5.55) are suppressed. So we simply set 
the metric fiuctuation at the boundary to zero Q\ 0. The solution B = ©o + ^ then follows 
from the matching procedure given in Section VB. We write, 

d{X) = iq^ X^) - G{X, X^)) , (5.58) 

where G'(X, X^) was constructed in (5.17). Plugging this answer into (5.54) we find at large 

4 = ^Qmdr {G{X,X^) - G{X,X^))\^^^^ ^ -^dr^GiX^,Xj , (5.59) 

where we have used the fact that at large separation G(X^,X^) is independent of r. Note 
that the coincidence limit G(X^,X^) gives the monopole self energy as discussed in Ap- 
pendix C. We have not carefully constructed G(X^,X^), however we can argue that the 
dominant radial dependence comes from the AdSs answer since the self energy is a local 
property which does not probe the full geometry. Reading the AdS^ answer from (C6), 
27rG(X^,X^) = lnln(f^/r^), we find the position of the monopoles in the saddle approxi- 
mation is, 

ln(f^/r^) qm 

Surprisingly this saddle exists even in the high temperature calculation. Examining the form 
factor R{k) (4.15) we see there are two terms which depend on in the exponential - the 
bulk to boundary propagator and the monopole self energy. These two terms work against 
each other, since the self energy wants to force the monopoles into the bulk. The radial 
integral is then dominated in a saddle approximation at the same as in (5.60). Note that 
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this contribution to the radial integrals is very different from the one we quoted in (4.18) 
which was dominated by a UV divergence. Roughly speaking in the Gn ^ limit there will 
be two contributions, one from the UV divergence given in (4.18) and one from the saddle 
at given in (5.60). Which one wins depends on the details of the UV completion. We did 
not see the UV divergent piece in the approach of this section since the saddle approximation 
does not account for contributions from the boundaries of field space (in this case oc). 
These do not minimize the action, but may still have the largest contribution. 

The precise meaning of this particular saddle value of is not clear to us. Physically one 
would expect the prefactor of the 2kF oscillations to be sensitive to the UV cutoff, however 
here we are claiming a universal contribution in the Gn ^ limit. 



VI. CONCLUSION 

We now briefiy summarize our results so far. We computed the contribution to the holo- 
graphic density-density correlator from monopole events in the bulk on a charged black hole 
background, dual to a field theory state with a charge density p. Essentially we computed 
the Witten-esque diagram in Figure 1 by placing a monopole-antimonopole pair in the bulk 
and integrating them with the appropriate action cost over the whole spacetime. Indepen- 
dent of the details of the bulk geometry, we saw that each monopole event had a Berry phase 
(3.9) associated with it arising from the coupling to the background electric field: 

where is the spatial coordinate of the monopole in the direction parallel to the horizon. 
This oscillatory phase directly translates into a singularity in the density-density correlator 
in momentum space at the precise wavevector = pq^j^. The precise form of the singularity 
depends on the monopole-antimonopole correlation in the bulk and thus on the geometry in 
question. We find at zero temperature: 

{5p{p)5p{-p))T=o - {vl{k - q^pf + ujlf+{k^ -k) (6.2) 
and at very high temperatures, 

mpm-p))T», - _ ^^py^ + ^ • (^-3) 

Here the unwritten prefactors of these expressions are both exponentially small and non- 
universal in that they depend on the UV cutoff of the theory. Vg is the speed of sound, and 
we have: 

K=—'^\ K'^ = —'^ (6.4) 
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A. The existence of a Fermi surface 



We now finally turn to an interpretation of these formulas. We begin with the location 
of the singularity, for which we must first understand the meaning of g^, the monopole 
magnetic charge. At the level of our current discussion it can be any magnetic charge 
that satisfies the Dirac quantization condition, qeqm = SttZ, where is the charge quantum 
associated with the compactness of the bulk U{1) gauge group. Note that this q^ has a direct 
boundary interpretation: it is the charge quantum appropriate to the global U{1) symmetry 
of the boundary theory Hilbert space. Now it is commonly thought [52] that in theories of 
quantum gravity such as AdS/CFT, the spectrum of allowed bulk charges should be such 
that all Dirac conditions are saturated^ i.e. the magnetic and electric charges should satisfy 

q^qm = 27r . (6.5) 

We will not discuss in detail the arguments leading to this (see e.g. [53]), but we assume it 
to be true. 

Then the location of the Friedel oscillation singularity in fc-space is 

fc. = ^ . (6.6) 

qe 

This is precisely the value of 2kF that we would find if we had constructed the entire charge 
density out of ordinary fermions with charge q^. Luttinger's theorem ^ interpreted here 
as a statement about a singularity in a gauge-invariant correlation function - is precisely 
satisfied^ despite the fact that there are no fermions in sight. There is something profound 
about the interplay between the field-theoretical Luttinger's theorem and the bulk Dirac 
quantization condition. It is also utterly clear that any deformation of the geometry will 
not change the location of the singularity: though p can be tuned, the relation between 
and p is topological. 

There are several ways to interpret this result. One is tempted to associate this singularity 
with the presence of a dual Fermi surface. Note however that we do not have (nor do we 
need) access to fundamental fermions in the gravitational description: thus what we have 
actually shown is not the presence of a Fermi surface so much as the presence of structure 
in momentum space at a special momentum that is related to the charge density 

by Luttinger's theorem. This conclusion, while not previously demonstrated in holography, 
is not surprising. There exist elegant nonperturbative formulations of Luttinger's theorem 
that essentially prove this statement [54, 55], although there to make the discussion precise 
one needs to explicitly use the UV structure of the Hilbert space, i.e. the degrees of freedom 
living on an underlying lattice. In those works one then further argues that if the system 
is in a Fermi liquid phase one should identify with the location of the Fermi surface. 

In a more general context it seems possible that the desire to identify this momentum with 
''the edge of occupied single-particle fermionic states" is not a useful one. In some sense 
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our holographic model provides us with a Fermi momentum with no reference needed to 
microscopic fermions. 

Of course in 1 + 1 dimensions such a construction is not entirely unexpected. In standard 
field-theoretical treatments there is ultimately little difference between a density of fermions 
and bosons in one dimension: independent of the statistics of the microscopic degrees of 
freedom, they result in the same low-energy effective theory (that of the Luttinger liquid), 
and they both exhibit Friedel oscillations at the appropriate wavevector. Our holographic 
model appears to be a different way to arrive at a similar model, although the low-energy 
theory is different, as we elaborate on below. 

We should note from general principle of gauge-gravity duality we expect the microscopic 
Hilbert space should have some large number ('W^") of gauge-charged fermionic degrees 
of freedom. If each of these fermions individually formed a Fermi surface, then one would 
expect to see Friedel oscillations at each individual wavevector, which would be of 0(1) 
in the large limit. An example of this was given in the recent field theory model in 
one dimension [56] of a "strange metal" - a gauge theory of Dirac fermions in the adjoint 
representation. The model has a single U{1) global symmetry, as in our setup, and a large- 
N limit can be taken. They find Friedel oscillations which preserve the original Fermi 
momentum of the underlying gauge-charged fermions, 2kjp ^ C^(l)- This is not what we 
find: our wavevector scales like the total charge density and thus is of 0{N'^). Somehow in 
the gravitational description any putative individual gauge-variant Fermi surfaces are not 
visible: however the existence of a single special gauge-invariant momentum that counts the 
total charge is essentially guaranteed by the arguments of [54], and reassuringly it appears 
in our holographic discussion. There is clearly more to be understood here. ^ 

B. Comparison to Luttinger liquid theory 

We turn now to more concrete considerations. In 1 + 1 dimensions the Luttinger liquid 
theory provides a robust field-theoretical framework for strongly interacting physics at finite 
density. Starting from a microscopic fermionic Hamiltonian, one linearizes the fermion 
dispersion relation near each Fermi point, obtaining a separate left and right-moving sector 
that form a 2d CFT. The theory can be rewritten in terms of the sound mode, treating it 
as a free fundamental boson 0, where interactions amongst the fermions modify the boson 
radius but are otherwise exactly marginal. The low-energy theory is thus characterized by 
two parameters: the speed of sound Vg and the boson radius, parametrized by K^ut- The 
density can be written in terms of this boson as (see e.g. [57]): 

PLutiy) = --V0(|/) + A (e^^'^^-e-^^^^^) + h.c. + •••). (6.7) 



One obvious difference with our model is at low energy and high densities they recover a CFT. We clearly 
do not have a CFT which would have separate left and right conserved currents. 
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The first term is the gradient of the sound mode; the second results in Friedel oscillations, 
where A is a UV cutoff indicating the non-universal nature of the term, and the • • • indicate 
higher harmonics. 

It is interesting to compare our results to those of Luttinger liquid theory. The form of the 
Friedel oscillations in a Luttinger liquid are easily calculated from (6.7). In the zero and high 
temperature limits they are of the same form as the our holographic results (6.2) and (6.3), 
with the substitution K K^ut^ "^s^Lut- In the Luttinger liquid these parameter 

relations are a consequence of 1 + 1 dimensional conformal invariance, which completely fixes 
the form of the finite-temperature correlation. However the low-energy physics of our system 
is not a CFT, due to the marginally running coupling k (2.24). This running of the coupling 
explicitly manifests itself in the expressions for K and (6.4), which contain extra scale- 
dependence relative to the Luttinger liquid case and can be thought of as measuring the 
value of at the scales corresponding to the density and the temperature respectively. We 
note that the precise dependence on the speed of sound is also different, again presumably 
due to its own logarithmic running (5.16). 

It is interesting to note that the essential role played by the monopoles in our discussion 
is to attach the second term (or some caricature of it) in (6.7) to the holographic current 
operator. Indeed in field theory that term may be thought of as the contribution of vortices 
(i.e. instantons in 2d) in the field to its two-point function. These vortices are very small, 
probing the UV structure of the theory, and so their effect on low-energy physics can be 
completely understood in terms of a pointlike local CFT operator e^^"^. In holography the 
story is somewhat different. The monopoles in 3d play roughly the same role as the vortices 
in 2d, except that they can move freely in the holographic direction, corresponding to a field- 
theoretical instanton with no well-defined size and preventing an interpretation in terms of 
a local CFT operator. Despite this fact, the dominant overlap of the monopole field with 
the high-momentum probe is when the monopole is very close to the boundary, i.e. when 
the dual instanton is as small as it can be: this results in a non-universal coefficient for the 
Friedel oscillations similar to that in (6.7). 

We note also that the zero-temperature compressibility of our model can be written in 
terms of K and Vg as (2.36) 

dp 
dfi 

where as above we have assumed that saturates the Dirac condition. This relation is 
precisely that of Luttinger liquid theory (see e.g. [57]), with only a single species of fermion. 
It is fascinating that this exact relation still holds in our large A^^ theory, except that we 
now have a K that is parametrically large, scaling with the number of charged degrees of 
freedom. In Luttinger liquid theory K can be related to the strength of interactions of the 
microscopic degrees of freedom; it would be helpful to understand precisely what this means 
in our model. 

We stress that while our model has some features in common with the Luttinger liquid. 



= (6.8) 

T=0 ^ 
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they are not equivalent: a Luttinger liquid is a pure CFT, and our marginally deformed 
theory must then be interpreted as describing a different phase of matter. Note for example 
that this system exhibits the familiarly distressing zero-temperature entropy associated with 
the near-horizon AdS2 x This entropy does not play an important role in our analysis. 
Indeed there appears to be very little relation between the Fermi momentum structure 
described here, which is intimately associated with bulk gauge field dynamics, and the 
recently constructed geometries thought to describe hidden Fermi surfaces due to their 
entanglement structure [58-60]. It would be useful to understand the connection (if any) 
between these two approaches. 

C. Future directions 

There are several directions for future work, some of which we highlight here as we tie 
up some loose ends in our discussion. 

1. Chern-Simons theory and UV completion 

As we have mentioned before, our usage of a Maxwell term in the bulk to describe a 
conserved current in two dimensions is somewhat nonstandard; as argued above this 
means that the theory is not conformal and will break down in the UV at a scale 
r\. One might instead be interested in studying a CFT, where a conserved current 
necessarily has separately conserved holomorphic and antiholomorphic parts. This is 
dual to the fact that it is represented in a 3d bulk by two gauge fields with Chern- 
Simons actions, i.e. 



In general L and R have Maxwell terms as well, but the dynamics is dominated by 
the Chern-Simons terms. Similar models were studied with an eye towards Luttinger 
liquids in [61, 62]. 

We note here that the dynamics of monopoles is very diflFerent in Chern-Simons theories 
[63]. A Chern-Simons term associates fiux with gauge charge, and thus a monopole 
creates charge as well as fiux. A monopole and an anti-monopole are then necessarily 
connected by the worldline of a charged particle, which can be viewed as creating a 
tension-full string that confines them [64]. So monopoles do not proliferate freely in 
Chern-Simons theories, and the physics discussed in this paper will not obtain in quite 
the same manner. 

We see that this construction is not obviously related to the single propagating Maxwell 
gauge field that we studied. However if one considers this bulk theory in a Higgs phase 
where a diagonal ("axial") combination of the gauge fields R — L is Higgsed, then it 
turns out that the resulting theory has a propagating 3d gapless mode that can be 




(6.9) 
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described with a Maxwell action for the orthogonal ("vector") current L + R (see e.g. 
[65]). This Higgsing can be accomplished by introducing a scalar field $ charged under 
the axial symmetry and explicitly sourcing it. The resulting physics is intricate, but 
if the scalar field fiows to a constant in the infrared then the IR physics will likely be 
described by the Maxwell theory studied here. This provides a possible UV completion 
for this theory, and the interplay between monopoles and the Chern-Simons term is 
likely to lead to interesting physics in the holographic context. 

2. Confinement in the bulk is dual to a charge gap on the boundary 

We turn now to a separate issue. In fiat space (with no background electric field) in 
(2 + 1) dimensions, the proliferation of monopole events in compact QED causes the 
gauge field to confine with a gap that is exponentially small in the monopole action 
[26]. It is natural to wonder what this means in the context of AdS/CFT: i.e. what is 
the dual of confinement of a gauge field in the bulk?^ 

In such a phase bulk electric fiux is confined into tight fiux tubes, each of which has a 
finite tension per unit length. Now a state with a U{1) charge in the field theory has 
electric fiux at the boundary; however in the confined phase the only way to arrange 
this is to force one of these fiux tubes to poke through the AdS boundary. The field- 
theoretical charge that we can get in this way is now quantized in units of the fiux per 
tube and each tube extends into the bulk, meaning that each quantum of charge 
costs a great deal of energy. The previous sentence is of course precisely the description 
of a field theory with a gap to charged excitations. Confinement of a gauge symmetry 
in the bulk is thus dual to the formation of an insulator in the boundary; a convenient 
representation of this theory would be the Sine-Gordon model in the bulk 



Note B has developed a potential, and fiux tubes in the bulk correspond to kinks in 

e. 

Now we should ask whether or not such a phenomenon happens in our model. Note 
that in our model the background electric field results in the familiar Berry phase 
that causes the Friedel oscillations. On length scales longer than the inverse density 
the contribution of a monopole-antimonopole pair to the functional integral oscillates 
wildly in their relative distance. Integration over their relative coordinate would appear 
to wipe out the amplitude, essentially leaving behind only "subtle" effects (such as 
the Friedel oscillations discussed here). 

However, one could imagine canceling these Berry phases by considering a ground state 
which spontaneously breaks translation invariance. This would correspond to a lattice 



Some of the ideas presented in this section were arrived at during discussions with R. Loganayagam. 




(6.10) 
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of flux tubes carrying the electric flux from the black hole horizon to the boundary, 
where it is interpreted as a discrete set of localized charges. However, the spontaneous 
breaking of translational symmetry in a (1+1) dimensional fleld theory is not allowed 
by the Coleman-Mermin- Wagner theorem; thus at flnite we do not expect this to 
be the ground state. More concretely, in the bulk there will always exist a gapless (in 
the fleld theory) Goldstone mode corresponding to fluctuations of this lattice, whose 
quantum effects are expected to destroy it. The resulting phase is what we have 
discussed in this paper: we do not believe that (6.10) is a useful starting point for 
understanding the flnite-density state. The situation is rather different if we imagine 
explicitly breaking translation symmetry with a lattice: here one should be able to 
understand conflnement (and hence holographic duals of insulators) at commensurate 
flllings. 

Note that most discussions of holographic "insulators" [66-68] describe a phase in 
which all excitations are gapped. Conflning the bulk gauge U (1) symmetry is different: 
only excitations that are charged under that U (1) is gapped, indicating that the charge 
dynamics are playing an important role. This is a very desirable property for the 
holographic description of an insulating phase. A neutral sector can remain gapless 
and mediate long-range forces between the charged excitations. It would be interesting 
to understand such phases further. 

3. Higher dimensions 

Finally, we turn to perhaps the most obvious extension of these results. As we have 
explained above, in 1 + 1 dimensions the appearance of a Fermi momentum without 
access to fundamental fermions is - while satisfying - perhaps not tremendously sur- 
prising: a similar phenomenon happens in the traditional fleld-theoretical treatment 
of one-dimensional interacting bosons. In higher dimensions, of course, there is a true 
distinction between fermions and bosons, and the subsequent appearance of a Fermi 
surface from a higher dimensional charged black hole would be correspondingly more 
interesting. 

The natural question then is whether one can include fluctuations of the appropriate 
magnetically charged excitations in higher dimensions to show that all charged black 
holes contain momentum-space structure. This is of course considerably more diflicult, 
as the dimensionality of the bulk excitation increases with dimension: in a 3 + 1 
dimensional bulk one needs to include monopole world-lines, i.e. a one-loop calculation 
involving quantum fluctuations of magnetically charged flelds. While this may be 
technically challenging it does not seem impossible, and we hope to return to this and 
the other issues raised here in the future. 

The recent paper [69] may provide useful guiding principles. In that work the electric 
flux through a bulk minimal surface is considered as an order parameter for fraction- 
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alization of charge. Indeed it seems that the observable we have constructed here 
is somewhat related to that order parameter. The monopoles play the role of the end 
points of the minimal surface (since the monopoles live in the UV of the geometry, this 
analogy is more precise than it sounds). The electric flux integrated along the surface 
gives the same answer as the Berry phase term - in our construction the path taken 
between the monopoles is not important because their are no explicit charges in the 
bulk (not behind the horizon.) In the language of [69] we are in a fully fractionalized 
phase. One big diflFerence is that the Berry phase here contributes an imaginary part 
to the action, in contrast to the conjectured real addition to the action of the minimal 
surface in [69]. It would be interesting to understand how the monopole story devel- 
oped in this paper changes with the addition of explicit bulk charges. We leave this 
to future work. 

We close on an optimistic note: the calculation outlined in this paper demonstrates 
that the class of charged black hole studied here behaves - when probed in the right way 
- as though it is made of charged excitations who are aware that quantum mechanics is 
encouraging them to form a Fermi surface. While there is clearly much more to be learned, 
we hope that the results presented here form another small step towards a useful application 
of holography to condensed matter physics. 
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Appendix A: Magnetic monopoles, electric charges, and Berry phases 

In this section we discuss some elementary aspects of the interaction of a magnetic 
monopole with a point electric charge in three spacetime dimensions. Some of our con- 
siderations can be found in a slightly different language in [70, 71]. In the bulk of the text 
we review the standard fact that the two actions 

Sa = J d'X^^F' , Se = J d'X^^iVef, (Al) 

describe equivalent theories, where the connection between the gauge field and scalar repre- 
sentations is (2.8): 

(rf^'^pg Q = l^F^P . (A2) 
9f 
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Sa is the natural choice for including electric charges, where the coupling to a charge moving 
along a worldline C is Sa ^ Sa + iqe Jc ^- Recall that this coupling tells us both how the 
fields respond to the charge (i.e. the Maxwell equations with sources), and how the charge 
responds to the field (i.e. the Aharanov-Bohm phase acquired along the trajectory C). Now 
consider working in fiat space in polar coordinates, 

ds^ = dr^ + dr^ + r^#^ (A3) 

with a point charge at the origin r = 0. Then the electric field is F^^ = i^^- Turning to 
the scalar representation we see from (A2) that changes as we wind around the origin 

ee(r,r,0) = ^. (A4) 

For B to be single- valued we should thus demand that it be a periodic variable ^ + q^] 
this will make sense for all possible charges only if the gauge theory is compact with charge 
quantum q^. 

Now consider adding a magnetic monopole source at a point X^, which in the Sq rep- 
resentation is the coupling Sq ^ Sq + iqmQ{Xjj^). If we now evaluate this coupling on the 
field Be produced by an electric charge we find the relevant term in the action to be 

Sherry — ^qm^i^m) — ~^ 0m • (^5) 

This phase associated with a magnetic charge in an electric field is precisely analogous to 
the usual Aharaonov-Bohm phase acquired by an electric charge in a magnetic field. Each 
monopole event should be weighted by such a phase in the functional integral. Note this 
action appears to depend explicitly on the value of the periodic spatial coordinate 0, but it 
should of course be single- valued modulo 27r; thus we conclude that we require 

qmqe = 27m n G Z . (A6) 

This is the Dirac quantization condition. 

The form of this Berry phase means that in some loose sense the wavefunction of the 
system after the monopole-mediated transition depends on 0, i.e. the monopole is creating 
angular momentum. We devote the rest of this Appendix to an understanding of this fact. 

First, we note that the canonical (Euclidean) stress tensor associated with Sq is not 
conserved in the presence of sources: we have instead 

v^T^" = gl (v,e [v^, V,] e + v^ev.e) . (a?) 

Here the first term is nonzero whenever B is multiple- valued, i.e. along the worldlines of 
electric charges. This is standard; demanding that this term (together with an appropriate 
contribution from the stress-energy of the particle itself) be conserved is equivalent to im- 
posing the equation of motion of the charged particle in a background field. However, the 
second term is nonzero on the monopoles themselves, where we find 

V^T^'^(X) = tq^\/MX)S^'\X - X^) = t5'^t>^5^'\X - X^) , (A8) 

ZTT 
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where the last equahty is for the single charge/monopole configuration studied above. As 
claimed, angular momentum is being created at the location of the monopole. Note the 
requirement that this change in angular momentum be integer is another route to the Dirac 
condition [72]. As described in Section VC2, this fact can lead to considerable confusion 
if one wants to couple the system to gravity, which is somewhat unhappy with a non- 
conserved stress tensor. The resolution that we adopt is that one cannot simply fix the 
monopole position and study the monopole-charge system in isolation: to perform a fully 
consistent treatment one must supply the extra momentum from somewhere, which in our 
problem is via a bulk-to-boundary propagator from the AdS boundary^^. 

We conclude by noting that this creation of angular momentum seems more transparent 
if we consider the U{1) gauge theory in a Higgs phase by condensing a scalar with charge 
— , i.e. the value of the charge that saturates that Dirac condition. In this phase the 
monopole operator simply creates a massive state, a superconducting vortex with fiux g^. 
It is well-known that the composite of an electric charge and a vortex has a shifted angular 
momentum [72] relative to the electric charge alone; thus the monopole operator must create 
it, just as the calculations above illustrate. 



Appendix B: Boundary term in variational principle 

Here we start with the bulk action in the form 

^2 



S=yI d'^VO i^Of + SMy (Bl) 

and construct Sbdy to achieve a variational principle that is consistent with the boundary 
condition (2.28). We note first that the on-shell variation of the bulk part of the action is 

SS,^ik = -gl [ d'yUSe, (B2) 
Jd 

However we would like the variation to be a function only of the source H{y) = U + ^V'^O. 
To this end we note that if we take 

SMy = glj^ dFy {^e - Y^i^Of^ (B3) 

then the on-sheU variation of the total action is 

SS = gl (fy e (^511 + ^ V^^^^ " / ^^^^ 



At some level the issues we are discussing - concerning a perfectly reasonable-looking source who never- 
theless violates a (contracted) Bianchi identity - are in many ways similar to the issues that arise when 
trying to deal with a magnetic monopole while taking very seriously the identity of the electric gauge field 
A. One might then speculate that there is an alternative way to formulate the 3d gravity problem that 
elegantly deals with it, or at least some analog of a "gravitational Dirac string" that one could attach to 
the monopole, but we will resist the temptation to explore this here. 
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as desired. Thus this is the correct form for Shdy Now the total contribution to the on-shell 
action (not just its variation) in the absence of any monopole sources can be computed to 
be 

So = fl^d'yeiy)Hiy) (B5) 

If bulk monopoles are added then the action will receive extra contributions localized in 
the interior. In the main text this contribution is evaluated on a monopole-anti-monopole 
configuration. In this case we use (3.18) for the field 9{x) as well as the relation between 
bulk-to-bulk and bulk-to-boundary propagators^^: 

G{r ta, yi] rs, 2/2) = K{r2, 2/2; Vi) (B6) 
to find the boundary contribution to the monopole-anti-monopole action: 

Se[X^, ^2] = IJ^ (fy H{y) (^q^{K{Xv^ y) - Ki^X^- y)) + j (fy'K{r = ta, y'; y)H{y') 

(B7) 



Appendix C: Computing monopole fugacities 

In this appendix we discuss the computation of the monopole fugacity (i.e. the action 
cost of a monopole as a function of position) in curved space. We first discuss precisely what 
information is captured by the fugacity. Consider a monopole at and an anti-monopole 
at Xjfi. It is shown in the text in (3.19) that the action of this configuration is formally 

^ ;;t^(2G(X^, X^) — G(X^,X^) — G(X^, X^)) (CI) 

Here the first term is an interaction energy, and the second two terms appear to be singular 
self-energies. These terms contain nontrivial physics: for example, the charge-neutrality 
of this configuration guarantees that any IR divergences present in the interaction energy 
cancels against a corresponding divergence in the self-energies. It is in this sense that 
we can write an IR-cutoff independent answer for the propagator in expressions such as 
(4.10). There is also UV information in these self-energies. In fiat (or presumably any 
maximally symmetric) space it is understood that the singular part of this self-energy is 
entirely accounted for by a constant term Sc that represents the action cost associated with 
the core of the monopole, and which can be computed given a UV completion of the theory. 

However if we break some of the spacetime symmetries then there is more information: 
while the UV divergence at the monopole location is the same (and thus there is still some 
action associated with the core) the interactions of the monopole field with the nontrivial 
geometry can result in a UV finite contribution that varies as we change the location X^ of 
the monopole. This information is contained in the function G(X, X^), and we require a 



This relation holds for any choice of boundary conditions. 

It can be helpful to imagine the energy of a point charge near an infinite conducting plate. 
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prescription to extract it, i.e. to perform the split: 

1^ lim G{X,X^) ^ -{S, + S,,if{X^)) = -S^{r) (C2) 

where Sint{Xjn) is the position-dependent contribution that we seek. This requires the 
subtraction of a divergence as X ^ X^, and it is important that this subtraction be 
unambiguous and introduce no extra spacetime dependence for the quantity Sint{Xjn) to 
have a meaning independent of the UV completion. 

To do this, we assume that we can incorporate the information regarding the core of the 
monopole in a function F(X^;X) = where X^) is the proper distance 

from to X and / is thus a function of a single variable, where we take /(^ ^ 0) ^ ^ to 
cancel the Coulomb divergence at the monopole core. The regulated interaction energy is 
then ^ 

SseifiXm)^-^ lim (C3) 

f{y) is not unique; there are different choices which correspond to different choices for the 
internal monopole structure. However it is clear that only the finite part of f{y 0) is 
important in determining Sint{Xjn)^ which is thus ambiguous only up to a constant that 
is independent of X^. Thus we have separated the spacetime dependence from the UV 
structure of the monopole. By requiring that the subtraction depend only on the proper 
distance we are essentially asserting that the internal constituents of the monopole, whatever 
they may be, obey the equivalence principle and have no interactions with other fields. 

We now implement this prescription for the case of a monopole living in pure AdSa. 
For pure AdSs the Green's function appropriate to the boundary conditions (2.22) can be 
calculated exactly in terms of Bessel functions to be 



rr' log 



(C4) 

Note that AdS^ is maximally symmetric, and thus the full radial dependence of the fugacity 
will come from the logarithmically running boundary conditions, i.e. the dependence of 
the propagator. We first estimate this radial dependence. 

We begin by noting that the portion of this propagator that contains the term in Kili 
does not contain any dependence, and so appears conformal. One can check that this 
portion actually corresponds to the full propagator corresponding to conformal boundary 
conditions 9{oo) = at the AdS boundary. It thus respects the full symmetry of AdS space, 
and so once it is Fourier transformed to position space it is a function only of the proper 
distance between the two points. Following our discussion above, we see that this part of 
the propagator corresponds to the function / defined in (C3), and if we simply drop this 
part and take r ^ we will obtain the self-energy. 
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The self-energy is then 



Sseifir) 



pdp- 



^t^9f Jo r2 log ' ^ 



(C5) 



We will perform this integral in a rather crude manner which nevertheless (see Figure 5) 

z^^ and at high 



captures the leading r-dependence. Note that at low p the integrand is / 

p\og[pr^ j 

p the integrand is cut off by the exponential suppression of the Bessel function at p ^ r. We 
thus simply integrate the low p asymptotic form from an IR cutoff at to p ^ r to find 



Sseif{r) - 



■/vioglog(^ 



(C6) 



In writing this down we have dropped an r-independent IR divergent term log log (^) coming 
from the lower endpoint of the integration; as discussed above these divergences will cancel 
if we are looking at a charge-neutral configuration and so are not included in the definition 
of the fugacity. This double logarithm is slightly more palatable if one thinks of it as the 
logarithm of the proper distance from to the monopole location. As e~^^^^f appears in 
the functional integral and we always have r <C r^, it is clear that the monopoles want to 
be at small r, i.e. towards the interior. 




FIG. 5: Numerical evaluation of monopole fugacity (C5) (solid line) compared with analytic ex- 
pression (C6) (dotted line), with = 200. As described above only the r-dependence can be 
computed within our framework, and thus an r-independent (and in fact IR-divergent) constant 
has been adjusted in the fit. 



Appendix D: Backreacted equations of motion 



Write the linearization of Einstein's and equations in the absence of sources as: 

1 



-R^. + Ig^. (r + 2- ^(ve)M + ^v^ev.e ) = o 



M 



(Dl) 
(D2) 
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where the superscript 1 means hnearize to first order. The three constraints are: 



2Ert = tcoEhl' - ikh",' + -^X. '^^ (2«fc^"r + i^Ehl) (D3) 

2/r 



2Er:c = pO' + ikK' - -iojeK' + ^ ikh\ (D4) 

^^rr - --hr " W^^^ + -T^h , —h, + — - ——6 (D5) 



The four dynamical equations are: 

M = + (/V + /) - f ^ + ^ '^K - - '^hl (D6) 



2£„ = -ff,^' _ /£!^ft'; + g,hl - '^0 (D7) 

2^.. = -r'K' - ^T^fK' - ^p'hl + ikp9 (D8) 
2^:cr = rV/^r'' + 3r//i"/ - p'/i", + ^pcj^^ (D9) 

Collecting the first order fields in one big vector = {9^ /i^, /i^, /i^, h^^' ^ h^J ^ h^^'} with 8 
entries we can write the above equations succinctly as a first order constrained system 

C^X' = {X'y = {HxYj X' (DIO) 

where ji = (r, r, x) runs over the three constraints (D3-D5). Finally in this language we can 
write the equation for the gauge invariant variables given in (5.5) as, 

where a = 1, 2. The dynamical equations for in (5.6) derives from the following identity 
for the matrix S (suppressing indices): 

= H^S - SHx + tC (D12) 

where 7^ is a 2 x 3 matrix and H^j is the 2x2 matrix given on the right hand side of (5.6). 
This equation and in particular the matrices i/^r and 7 can be derived by brute force. Note 
that 7 is never needed so we do not write it here. 



Appendix E: Hydrodynamic limit 

We started studying the hydro limit in Section VA2 where a sound mode was found 
(5.14). Here we would like to go a little further and derive the width S of the sound mode. 
We will only consider T = here, although it is certainly possible to include a small T. 
Firstly we will need to go beyond the normalizable solution (5.13) of the hydro equation 
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(5.12) . So to begin with we find the second solution to (5.12). This can be done using 
standard tricks. The Wronskian of two solutions to (5.12), (cr^,cr^) and (crf,cr^) has the 
following form 

W = {a^at - a^a^) = WoD,{r) (El) 

where Wq is a constant and Di{r) = k'^f^r + 2r^CL;|;. We can then use (El) to integrate the 
unknown solution. The general solution is 



{a + bt{r)){^ \^,,.A^h-^\\ (E2) 





where a, 6 are integration constants and t{r) is defined as, 

= t{r ^ r,) = (^^) ' + ... + X (r - r,)° + . . . (E3) 

where the last condition fixes the integration constant in t. Once we have a, b we can simply 
read oS the boundary quantity /3/a which goes into the final form of the current current 
correlator (5.19) 

/?/«= .fi^rf', . (E4) 

^ cj|(a/6 + t(oc)) ^ ^ 

We now set out to find b/a. In order to find a width we need to effectively allow the 
sound mode to "fall into the dissipative horizon." This happens at any non-zero uje since 
the hydro approximation to (5.6) necessarily breaks down at (r — Th) ^ \(^e\' We should 
then match onto another solution close to the horizon. This procedure was used heavily in 
[17] which we follow here closely. The AdS2 x K region makes this matching particularly nice 
and has an interpretation in terms of RG fiow of double trace deformations in the boundary 
field theory [73-75]. 

The matching procedure requires splitting the radial direction up into inner < (r— r/^) <C 
Th (close to the Horizon or AdS2 x K geometry) and outer regions (r — Th) ^ \(jOe\ with 
solution given in (E2). There is an overlap region where we match if \(jJe\ <C Vh- We also 
take (jJe ^ k. It is easiest to proceed using second order equations for which we can be 
derived from the first order equations (5.6). In the inner region this equation becomes: 

-'/(C) + ^(^^-i(C) - ^-i(C) = (E5) 
where C = (^ ~ '^h)/^E' The regular solution at the horizon r ^ is, 

ai(C) = e-k'^<--\l + 2sgnK)C + 4C') (E6) 
We can now expand ai at large C and expand (E2) at small r ^ and match. 

-i(C) = 4C' + ^ - ^C-'sgnK) + 0{C') (E7) 

ai(r) = 2a(r - mf (1 + 0{r - r,)) - . (1 + 0{r - r,)) (E8) 

3 (r - Th) 
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We can clearly match the two power laws (r — ThY and (r — rh)~^ (which incidentally fix 
the conformal dimension of the appropriate operator in the AdS2 theory.) There are various 
terms in the above expansions that do not match, these were understood in [17] as coming 
from perturbative corrections to both the inner and outer solutions. They are generally 
associated with certain analytic corrections to the real part of the sound mode dispersion 
relation, and will not effect the width [17]. So we can ignore these and accurately read of 
the ratio, 

6/a= 4KA/^)'sgnK) (E9) 
lo 

Putting everything together (E9), (E4) and (5.19) we find the width as defined in (5.14) to 
be, 

'8rl \n{n/rh) 

Continuing this result to real frequencies we take only the ue > section and continue S(c<j) 
in the upper half plane to get the retarded function 



E = a;^sgnK)A_^^^-^ (ElO) 



^{ue ^ -lUj) =^R = ^ / \ (Ell) 

8r^ln(r^/r/,) 

For > Th it is clear that the sound mode appears in the lower half complex frequency 
plane as required by stability. 
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